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Abstract 



This thesis is a study of two dimensional noncritical string theory. The main 
tool which is used, is the matrix model. 

There are several chapters. After a general introduction there follows 
an introduction to the Liouville model where the fundamental issues of its 
formulation are discussed. In particular, the special states are introduced. 
Then, in chapter three, some calculations of partition functions on genus one 
are given. These use field theory techniques. The results are compared with 
the matrix model. In chapter four the matrix model itself is introduced. 
Some of the concepts and relations which are used in later chapters are ex- 
plained. Chapters five and six include comments on two important subjects: 
nonperturbative issues and string theory at finite radius. Chapter seven is 
devoted to zero momentum correlation functions as calculated in the matrix 
model. One important result is a set of recursion relations. Chapter eight 
extends the treatment to nonzero momentum. The main result is a clear iden- 
tification of the special states. The chapter also includes some comments on 
the Wheeler de Witt equation. Chapter nine introduces the matrix model 
Woo algebra. This organizes the results of previous chapters. In particular, 
a simple derivation of the genus zero tachyon correlation functions is given. 
Chapter ten extends the results of chapter nine to higher genus. It is seen 
how a deformation of the algebra is responsible for much of the higher genus 
structure. Some very explicit formulae are derived. Then, in chapter eleven, 
the Liouville and matrix model calculations are compared. Finally, chapter 
twelve is devoted to some general conclusions. 

The advisor for this work has been Prof. David Gross. 
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Chapter 1 
Introduction 



The outstanding problem of todays physics is the unification of quantum 
mechanics and general relativity, quantum gravity. It is widely believed that 
such a theory would also provide a unification of all forces and therefore give 
a unified picture of all physics. The putative new theory would hopefully 
also predict a multitude of new and fascinating phenomena. Needless to say, 
such a theory does not exist today. 

During the past decade or two however, a very promising candidate, in 
fact the only existing candidate, has been developed. This is string theory. 
In string theory several of the problems inherent in quantum gravity seem 
to be resolved. In particular the nonrenormalizability of the perturbative 
expansion is taken care of by a Planck scale cutoff given by the size of the 
string. 

Unfortunately a phenomenologically acceptable string theory does not 
so far exist. The heterotic string, the most promising candidate, gives a 
reasonable gauge group which by symmetry breaking can be reduced to the 
low energy SU(3) x SU{2) x £7(1) which we see. However, there does not 
seem to be any unique way of compactifying the 10 space time dimensions 
typical of a critical super string, down to the four dimensions of the real 
world. For a general introduction, see [0, and references therein. 

If we believe in the promise of string theory, but are discouraged by the 
lack of immediate success of the critical theories, we should consider alter- 
natives. Such an alternative is provided by the noncritical string theories. 
These string theories rely on solving two dimensional quantum gravity which 
in some cases turns out to be a feasible exercise. This invalidates the reason 
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for restricting oneself to the critical dimensions. Unfortunately, technical 
obstacles stop us from directly constructing a four dimensional string the- 
ory. The best we can do so far, is to obtain a string theory describing a 
string moving in a two dimensional space time. This theory lacks a lot of 
the complications of higher dimensions but is complicated enough to contain 
interesting physics. 

Remarkably there also exists a new tool, the matrix models, which allow 
us to solve these theories exactly to all orders of string perturbation theory. 
This two dimensional string theory, in particular its matrix model version, 
will be the subject of this thesis. 

Given the simplicity of its formulation, it is very surprising that the two 
dimensional noncritical string theory was not studied long ago. While the 
precise connection between this theory and the Liouville mode quantum grav- 
ity approach requires some technical machinery, there are no difficulties es- 
tablishing it as a perfectly consistent two dimensional string theory. 

The important point to realize is that a nontrivial background, in this case 
a linear dilaton field, can change the effective central charge and hence the 
critical dimension. This is of course just the string theory way of describing 
the conformal field theory construction with a background charge. 

More generally from the string theoretic point of view, one needs to solve 
the j3 equations. These are precisely the requirement that the world sheet 
theory is scale invariant, no conformal anomaly. Without a dilaton field they 
reproduce the Einstein equations for the metric. This is true at tree level for 
the world sheet a model. In general, there are corrections at higher orders 
in a', where \fa' ~ size of the string. In two dimensions we really need the 
nontrivial dilaton background to get a consistent theory. But as was shown 
in [55], ^] this is not the only possibility. There is a one parameter family of 



solutions describing a black hole background, the parameter being the mass 
of the black hole. 

Clearly the two dimensional string theory can serve as a laboratory for 
testing many of the ideas of string theory. It might give some insight into 
the deeper issues of how string theory gets around the problems of quantum 
gravity. In the end this could turn out to be the single most important 
contribution of string theory. Even if string theory turns out not to be the 
correct theory of quantum gravity realized in nature, we might learn a lot by 
studying how it resolves the quantum gravity enigma. If nothing else, this 
solution, even if it is wrong, can serve as an important source of inspiration 
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in the work towards the correct theory. 

There are many mysterious aspects of quantum gravity and it is so far 
unclear in what sense they are resolved by string theory. One of the most 
interesting ones is the problem of black holes and the loss of quantum co- 
herence. It appears as if important quantum phase information can get lost 
without a trace into a black hole. After all, it is claimed that a black hole has 



no "hair" ||12|| . As a result a pure state can evolve into a mixed state. This 
breaks the unitary time evolution of quantum mechanics. It is important 
to note that one cannot simply say that the information is there somewhere 
inside of the black hole. Eventually the black hole may evaporate due to the 
Hawking process revealing nothing. The Hawking radiation is also supposed 
to be purely random and thermal containing no information. How is this 
mystery solved? 

It has been proposed that black holes in fact do have a lot of quantum 



hair [12]. This is not forbidden by the classical no hair theorems. Clearly 
this would have the potential of solving the quantum coherence problem. It 
has also been proposed, ||29|| , that there is extra hair in the two dimensional 
string theory which we are about to study. The hair would be associated with 
the Woo symmetry which we will encounter at several places throughout this 
thesis. If this is true also for higher dimensional string theories as claimed in 
p9[| , string theory would have solved one more of the fundamental issues of 
quantum gravity. In some sense this would be an even more remarkable and 
far reaching achievement than just taking care of the nonrenormalizability. 

However, so far none of these speculations have any secure basis. In fact, it 
may very well be that the solutions to problems like quantum coherence have 
nothing or very little to do with strings. In particular it is conceivable that the 
loss of quantum coherence is just what it is. A sign that quantum mechanics 
as we know it breaks down, as suggested by Hawking Depending on 

taste, this could be an even more exciting perspective. Before getting to 
deep into string theory, let us briefly remind ourselves what this could mean. 

We have come to accept the strange duality which exists in quantum me- 
chanics as a rather fundamental aspect of nature. The unitary time evolution 
on the one hand, and the somewhat magical, certainly nonunitary, collapse 
of the wave function at the moment of observation on the other hand. The 
second process defies a well defined physical description. When does it really 
take place? The philosophical discussions about this seem endless. It would 
clearly be desirable to have a physical understanding of this. Is it then too 
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farfetched to suggest that a nonunitary quantum mechanics, the nonunitar- 
ity supplied by gravity as suggested by the loss of quantum coherence above, 
could give a physical process for the collapse of the wave function? Rather 
than just trying to device ways of getting rid of the apparent nonunitarity, 
we should carefully ask ourselves whether we might not need it after all. The 
point to be made is that it is very likely that string theory may provide a 
solution to the problem of reconciling quantum mechanics and gravity, but 
the question is whether that is what we want to do! Is this really the solution 
which nature has chosen? 

Clearly the quest for quantum gravity is a very open one. Even if string 
theory is our best candidate for the moment, we should be open minded. It 
is not obviously true that strings really have the capacity to solve all deep 
conceptual problems which confronts us in e.g. quantum gravity. However, 
string theory clearly deserves the attention it receives for whatever clues it 
might give. With this in mind let us embark on the careful study of two 
dimensional string theory which is the subject of this thesis. 
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Chapter 2 

The Liouville Model 



2.1 Introduction 

In this chapter we will review the two dimensional noncritical string from 
the field theoretic point of view. In a later chapter an introduction to the 
matrix model approach will follow. 

We will also consider the theory as a model of two dimensional world sheet 
quantum gravity. For this we will briefly review some different approaches. 
Historically this was also the way in which one first solved the model. Only 
later the focus has been more on the string theoretic interpretation. Ironically 
it is really only the quantum gravity interpretation which need a nontrivial 
justification, as we will see. 

Towards the end of the chapter we will discuss the "special states" . These 
are probably one of the most interesting aspects of the theory. Since they 
are remnants of the usual excited string modes in higher dimensional theo- 
ries, they could potentially tell us a lot about stringy phenomena. We will 
encounter them again and again throughout this thesis. 

2.2 The Way to Solve it 

Our focus will be on the gravity coupled c = 1 model and its target space 
interpretation as a two dimensional string theory. We will however begin by 
considering the more general case of a minimal model with central charge 
c < 1 coupled to gravity. More precisely we will be discussing induced 
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Liouville gravity. 

The action for a minimal model coupled to Liouville gravity, in complex 
coordinates z = o\ + io-i and conformal gauge h = p ^ ^ ^ ' * s ^* ven ^ 

J {dXdX + l -^VhRX + <90<90 - Q-VhR<p). (2.1) 

X is some matter field, while <ft, as we will see in a moment, is the grav- 
itational field. Recall that in two dimensions the metric has only one in- 
dependent component. R is the world sheet curvature while a and Q are 
background charges for matter and gravity respectively. Due to the "i" in 
the action, the model is unitary only for some special values of a. Good 
introductions to these and other related issues in conformal field theory may 



be found in |TT] and [|31|]. a is adjusted to give the appropriate value of the 
central charge of the minimal model we are considering, while Q is tuned in 
such a way that the total central charge of both matter and gravity adds up 
to 26. Hence it is also possible to interpret this as a consistent string the- 
ory. X and <p are then some target space coordinates. The curvature term 
describes the coupling to a linear background dilaton field. 

Given the action Q2.1|) the stress energy tensor, defined by T af3 = -^Tjff^, 
is easily seen to be 

T zz = --(: dXdX : -iad 2 X) - -(: d<pd<p : +Qd 2 <p) (2.2) 

by varying the action with respect to the metric. The terms linear in the 
fields come from varying the world sheet curvature and some following partial 
integrations. From this the Fourier modes of T zz , i.e. the Virasoro generators 
LnS, are extracted as 

1 00 1 i 

L n = o /J : ( a n-m<y m + Pn-mPm) ■ -^ a { n + l ) a n ~ ^Q( U + l )Pn- (2.3) 

ot n and p n are the matter and gravity oscillators respectively. In particular 
we have 

L = \{P 2 + V\ ~ ®P ~ iQp<t>) + ( a -rn®m + (3-mPm)- ( 2 -4) 
^ m>0 
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The energy momentum tensor can also give us the relation between a, 
Q and c, which we need. For simplicity of notation let us just consider the 
matter part. The propagator is given by 

< X(z)X(w) >= -\og(z- w), (2.5) 

and the corresponding expression for the antiholomorphic part. It is a simple 
exercise using Wick contractions to derive the operator product 

r(W ^f"? +- (2.6) 

from which one can read 

c = 1 - 3a 2 . (2.7) 
Similarly one finds for the Liouville part 

c L = l + 3Q 2 . (2.8) 

Since we need c + Ci = 26, this fixes Q to 



Q = \ 2± r^- (2-9) 



In the original works, |[T], ^7, [52], ftjfl , this action was mostly interpreted 
as describing two dimensional world sheet quantum gravity coupled to some 
matter. Only later it became fashionable to think of it as a noncritical string 
theory. The latter interpretation is clearly justified without further calcula- 
tions. We have achieved quantum scale invariance, the hallmark of a string 
theory, and later we will also see how to construct scattering amplitudes. To 
really make the quantum gravity connection, some further work is however 
needed. Let us briefly review the different steps. 

The starting point is to note that the classical scale, or Weyl, invariance is 
broken quantum mechanically. The need to regularize introduces a potential 
dependence on some scale which only goes away in the critical dimension. 
The breaking of scale invariance is measured by the trace of the energy mo- 
mentum tensor, T zS . Let us for simplicity focus on c = 1 where there is no 
background charge. Classically one would certainly expect the trace of the 
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energy momentum tensor to be zero, but the quantum mechanical answer is 
in fact 

T " = ~Ya pR (2 - 10) 

where R is the world sheet curvature. Let us give a simple, heuristic but 
perhaps illuminating derivation of this result. 

The key point is the properties under coordinate transformations of T = 
T zz (or T = T zz ). For c = 1 the matter part of T zz is just — \ : dX(z)dX{z) : 
Due to the normal ordering a coordinate transformation z — > z(z) has the 
following effect: 

--lim(d z z{z + e/2)d z z{z-e/2)dX{z{z + e/2))dX{z{z-e/2))-\). (2.11) 

The e limit and subtraction is the conventional point splitting version of 
normal ordering. Some simple algebra then leads to 

T(z) = (d z ~zyf(z(z))-±{~z,z}. (2.12) 

The last term is the Schwarzian derivative defined by 

r „ , d 3 z 3 fd 2 z\ 2 /010 , 

We now try to reproduce this transformation property classically by adding 
an extra term to T zz = —jdXdX. This term must be 

-^{/.*> ( 2 - 14 ) 
for some function / = f(z,z). It is easy to show that indeed 

{f,z} = (f/{f,z} + {z,z}. (2.15) 

We next need the conventional law of conservation of the energy momen- 
tum tensor 

V- z T + V z T z - z = (2.16) 
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The extra term in T must be accompanied by a similar extra term in T zl . 
Otherwise (|2.16|) would in general hold not hold. It is easy to see that an 
extra term in T zz of the form 



-^ddlogdf (2.17) 

would do the job. But this is just — jgpR in conformal gauge where R = 
— 2p~ 1 dd\og p with the conformal factor p = df\ 

In the case of general c the expressions are multiplied by the central charge 
c. This can be seen from (|2.6|) which upon contour integration against the 



parameter / produces the infinitesimal version /'" of the Schwarzian deriva- 
tive. The net effect is that we have generated gravitational contributions to 
the energy momentum tensors given by 

c 



T^ z = --{f,z] (2.18) 



and 



Tg = ~pH (2.19) 

That R is the only thing that T zz can be proportional to is really clear 
already from dimensional analysis and coordinate invariance. The precise 
coefficient, however, is provided by the above argument. We will come back 
briefly to the different expressions for the energy momentum tensor above 
when discussing some different approaches to two dimensional quantum grav- 
ity, but let us first see what to do with fl2.10|) . 



The central charge in ( 2.1U| ) receive contributions both from the matter 



fields and the ghost fields. This gives 

26 — r 

T z - Z = A^PR (2-20) 

To get scale invariance without c = 26 we in some sense need R — 0. The 
consistent way to achieve this is by thinking of R = as an equation of 
motion for the two dimensional metric. As an aside one may note that 
this "induced gravity" is different from ordinary Einstein gravity. There the 
action is simply J \f~hR which in two dimensions is a topological invariant, 
the Euler characteristic. The action of the induced gravity is instead, in a 
general gauge, 

26 — c r i i 1 

J d 2 xd 2 y^h(x)^h(y)R(x)-R(y). (2.21) 
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Hence nonlocal! After some initial despair one finds that the action is local 
in conformal gauge which therefore, clearly, is a very sensible choice. We 
then get 

26 — c r /T 1 

— — / dVYMo^W + m (2-22) 

where p = with the metric on the form h = ph. h is some fixed background 
metric, e.g. with constant curvature. For convenience, one usually rescales 
the field <fr to obtain a more standard kinetic term 

i- J d 2 a\fi(d a (j)d a( f> - Q C R<P) (2.23) 

where 



Qc = \j\^- (2-24) 

Henceforth we will drop the hats. In string theory language we have fixed 
the a' for the Liouville coordinate as a' = 2. So far so good. Now the tricky 
part comes, which is the reason for why the model was not solved until quite 
recently. The action above is the classical action. To quantize we need to 
do the path integral. A path integral needs a measure, and the measure 
needs the world sheet metric for its definition. But the world sheet metric 
is a dynamical variable! Clearly a very confusing situation. The trick is to 
make a change of variables in the measure and define it with respect to a 
new fixed background metric. The change of variables presumably leads to a 
Jacobian. The assumption in [27] was that the only thing that happens is a 
renormalization (finite) of the different couplings in (|2.23|) . In particular, one 
needs, for consistency, that Q c as given by fl2.24|) renormalizes into Q as given 



by ( |2.9D . This change of 26 — > 25 can also be understood as coming from the 
new quantum contribution to the central charge from the Liouville mode. 
Independence of the background metric translates into the requirement of 
scale invariance in the theory ( |2.1| ). 

From many points of view the argument above is purely hand waving, 
although it is hard to imagine what else could happen. More rigorous ap- 
proaches however do exist and are described in the literature [^3] . 



Although this is the common way to do things nowadays, the first solu- 
tion was not formulated in the conformal gauge but in a peculiar left-right 



asymmetric gauge in |52], |>9| . For completeness we will give a short descrip- 



tion of this. This will also throw some further light on the expressions for the 
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energy momentum tensor. The treatment will be very brief and interested 
readers are referred to [52, 59 1. The reasons for including this section are 
mainly historical. 



2.3 Another Solution 

As we already have seen the energy momentum tensor of the Liouville theory 
is given by 



T 



, ; 26- c(d 2 p 3 fdp^ 2 



12 

T G + = -^-«9_<9 + logp. (2.25) 




Both these expressions can be obtained by varying (|2.22j ). We have temporar- 
ily changed to a Minkowski metric to agree with the conventions in [[52], ^ . 

If we write p(x~,x + ) = d^f(x~,x + ), the T component is the Schwarzian 

derivative of a new field /. Some simple studies of the coordinate transforma- 
tions involved show that / can be thought of as a coordinate transformation 
which connects conformal gauge with a gauge where the metric looks like 

(2.26) 

\ '"++ j 
with 

h ++ {f{x-,x + ) 1 x + ) = -d + f (2.27) 

This is the left-right asymmetric gauge used in |]52| , |69|| where the Liouville 
theory was solved for the first time. 

A crucial ingredient in this solution is an apparent SL(2, R) symmetry. 

One way to see this symmetry is to note that T is invariant under precisely 

such coordinate transformations in /, with coefficients being arbitrary func- 
tions of x + , since it is a Schwarzian derivative. In other words, invariance 
under 

^ a{x+)f + b{x+) 

J c(x+)f + d(x+y { ' J 

The SL(2, R) structure gives rise to Kac-Moody currents which are help- 
ful when one attempts to solve the theory. They are the following components 
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of the metric field h ++ 



h ++ (x~,x + ) = J+(x+) -2J°(x+)aT + J-(x + )(x-f (2.29) 

subject to the equations of motion. Note that the equations of motion, 

T = 0, imply, but are not equivalent to, dj-h ++ (f,x + ) = 0. One can also 

show that the components J transform as they should under the symmetry. 

The weights of the SL(2, R) algebra is then related, through Ward iden- 
tities, to the gravitational dimensions. The results obtained agree with the 
conformal gauge approach. 

One way to understand the difference between the two different gauges 
is to note the difference between constraints and equations of motion. In 
conformal gauge, the equation of motion is T_ + = 0. In the left-right asym- 
metric gauge it is instead the less restrictive T =0. T_ + = is in this 

case a constraint which we finally have to apply to get the correct physical 
states. They are hence of highest weight with respect to the SL(2, R). 



The basic trick of [[52], is not so much the different gauge choice. The 



SX(2,R) symmetry can be faked directly in the conformal gauge [|13[]. One 
only needs to make a field redefinition of the form p = d-f. Substitutions 
like this are something one can contemplate in more general theories. The 
resulting theory has higher derivatives and looks very different from the orig- 
inal one. Equivalence is guaranteed in the end only by applying appropriate 
constraints. Two dimensional gravity is an example where the higher deriva- 
tive theory in some sense is simpler to study, thanks to the new symmetry 
appearing. 

In the following we will leave this method of solution aside and stick with 
the conformal gauge treatment which has become the standard way to do 
things. 



2.4 Correlation Functions 

Let us now limit ourselves to the case of main interest in this thesis, c = 1. 
From formulae (p.7| , |2.9|) we see that a is just zero (obviously) and Q = 2y2 
in the usual convention with a'^ = 2. For reference we write down the action 
again 

— J (faVh(—d a Xd a X + ^d a <pd a <j) - V2R4>). (2.30) 
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Let us now see how to construct the operators of the theory. As usual in 
string theory, states are created by vertex operators sitting on the surface. 
These need to be integrated to yield string scattering amplitudes. Con- 
sistency therefore requires that they have conformal dimension (1,1). We 
therefore need the L generator in (|2.4|) for c = 1. 

If we apply this to a state without oscillators (the general case will be 
studied later) i.e. a tachyon, we obtain the condition 

^f-^(2V2 + k) = l (2.31) 

with solutions 

k = -V2± ^l^ 1 . (2.32) 
v2 

Compared to ( |2.4|) we have redefined p^ = —ik. We see that there are two 
solutions. Quite surprisingly these two different solutions are of very different 
character JF4J . We will come back to this later. For now we simply choose the 
positive sign and denote this choice as the "right dressing". The "wrongly 
dressed" states will be briefly discussed elsewhere. 

The explicit calculations of tachyon correlation functions are really very 
simple . Essentially, one may borrow the well established results from the 
critical string. Things are in fact even simpler here, and a lot more explicit 
calculations can be done. 

The main new ingredient, which also is responsible for many of the pecu- 
liar features of two dimensional string theory, is the very restrictive kinemat- 
ics. Let us assume that we want to compute an N-point tachyon correlation 
function. To do so, we need to impose momentum conservation. The matter 
part is just the standard 

N 

5> = 0. (2.33) 

i=i 

The Liouville part receives a correction due to the coupling to the world sheet 
curvature and becomes 

N a/t7 I r>- I 

£(-V2 + ^r 1 ) = 2V2(g - 1). (2.34) 
i=i V2 

If we furthermore assume that 

Pi < 0, and pi > for i ^ 1 (2.35) 
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we find 

'a' | pi 



-V2(0-l + N/2). (2.36) 

The momentum of the single negative chirality tachyon is hence fixed! 

An important complication arises if there is a world sheet cosmological 
constant. In string theory we would say that there is a nontrivial tachyon 
background which the string couples to. The form of this term is determined 
by conformal invariance to be 

Ae^ (2.37) 



where a = — a/2 for c = 1 and a = — ^ + y — 2 for a general minimal 
model. This follows from ( 2.31|) , the condition that the operator is (1,1). 



This means that Liouville theory with a cosmological constant is not a 
free theory on the world sheet, but in fact interacting with vertices of infinite 
order. Even though a straightforward solution of the theory which produces 
all amplitudes unambiguously has ,as yet, not been obtained, there exist sim- 
ple arguments and analytical continuations which give all amplitudes. The 
best verification that these arguments make sense is of course the agreement 
with the matrix model. 

The basic effect of the cosmological term is that it can modify the mo- 
mentum conservation law ( [2.34 ) above. One insertion of the cosmological 



term injects an amount a of Liouville momenta. Hence, if we break ( |2.34| ) 
by q the amplitude can still be nonzero and will be proportional to 

A q/a . (2.38) 

This makes perfect sense in the case of n = q/a an integer. Provided we 
by hand insert n extra punctures we can make all computations in the free 
(and X) theory. We must of course remember also to do the zero mode 
integration 

d a a 

which give the dependence on A. The left hand side is obtained using 
J JhR = 8tt(1 - g). We note that for genus one, without any insertions, 
Liouville momentum is conserved and the integration gives rise to a volume 

I — log A I . (2.40) 

a 
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This will be used in the next chapter. 

A serious problem, however, arises for noninteger n. How do we insert 
a fractional number of punctures? Luckily it is possible to argue, as has 



been done in p3[ , that the full answer can be obtained by "analytically" 
continuing the integer answer. 

Although we should be grateful, and perhaps surprised, that it is possible 
to get away with a trick like this, it would clearly be more satisfying if 
the calculation could be done directly. This, however, would require an 
understanding of path integrals beyond the present state of the art. Unless, 
of course, one uses the matrix model. 

For future reference, let us give the explicit expression for the tachyon 
correlation function as obtained in []23|, [5B|. In principle the calculation is 
very simple. Just compute the Veneziano like integral with the N tachyons 
and the appropriate number of screening charges. The answer is 



nffi^o^^-. (2 . 41) 

l h rW|p|) V" 3 



The general integral is in fact a bit tricky and one needs some analytical 
argument to easily derive it. 

In the above equation we have changed A to its Legendre conjugate \i. 
This means that we are considering 1PI rather than connected amplitudes 
with respect to the puncture. In chapter 4 we will consider this in more 
detail and find, among other things, that /i ~ F° r generic momentum 

the transition is hence very simple. One simpfy absorbs the external leg 
Y^-^ for each screener into A and replaces it by /i. The is provided 
by the p^y associated with the zero momentum punctures in analogy with 
(|2.41|) . At zero momentum we should also remember to amputate external 



zero momentum tachyon legs. For, e.g., the two point function we then get 
logyU rather than j^j- We will come back to this on several occasions where 
we also will see that one must be much more careful when at zero momentum. 



2.5 The Special States 

If we regard the c = 1 quantum gravity as a critical string theory in a two di- 



mensional space time with the Liouville field as the extra dimension [39, 55], 
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the naive expectation would be that the massless tachyon completely ex- 
hausts the spectrum. A simple light cone gauge argument would indicate 
that, since there are no transverse dimensions, there are no physical excita- 
tions except for the center of mass of the string-the "tachyon". Indeed, in 
the computation of the one loop partition function (i.e. genus 1) using con- 
tinuum methods as in the next chapter, only the tachyon seems to contribute 



[y, (14j| . However it turns out that for certain discrete values of the momenta 
there exists new states and new nontrivial operators. In the context of the 
matrix model these states were first seen in the calculation of the puncture 



operator two point function PD[. They also appear in the external legs of 



multi point puncture operators |3(| |69fl . In conformal field theory it is well 
known that such special states appear for c = 1. In this section we review 
this story and show how the special states appear in c = 1 conformal field 
theory. 

From the point of view of string theory it is not surprising that there are 
other degrees of freedom. After all, if this theory is a two dimensional theory 
of space time the two dimensional metric should be a dynamical degree of 
freedom. In two dimensions there are of course no propagating gravitons, yet 
there are global, or topological, degrees of freedom associated with the metric. 
Also, in string theory the string coupling is a dynamical degree of freedom. 
This corresponds to the zero momentum component of the dilaton. The 
full set of special states, with nonvanishing momentum as well, presumably 
correspond to topological degrees of freedom of the two dimensional string 
theory. They are the physical remnants of all the massive modes of the string 
in higher dimensions. If we wish to be able to construct the most general 
solution of two dimensional string theory we must be able to excite these 
modes. 

A c = 1 conformal field theory has extra primary states which are closely 
related to the null states of this theory. Null states are zero norm states that 
are created by descendents of the primary field e tpX for special values of the 
momenta p. These special values are quantized in units of The relevant 
primary fields are labeled by two integers, r and s, as in the minimal models. 
The momenta take the values p = ?-j=, and the field e ipX has conformal 
dimension 

= (r-sl = a V (2A2) 
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The null states are descendents at level rs, hence have conformal dimension 



(r + s) 2 



''null 



(2.43) 



The states have both a left and a right moving component, with h = h, but 
we will usually suppress the left moving component. In the minimal models 
all primaries have descending null states. This is not the case for c = 1 where 
the momenta, p, in the uncompactified case, need not be equal to the values 
given by fl2.42|) , but can take any value. 



The reason for the existence of new primary states is that the null states 
are not only null (zero norm), but in fact vanish identically. The map be- 
tween states and conformal transformations, i.e. the Virasoro generators 
L_ n , degenerates. Hence there exist new states which cannot be obtained 
by conformally transforming the primaries above and are therefore new pri- 
mary states. Recall that the descendants of a primary are precisely those 
states which can be generated by acting with conformal transformations on 
the primary. One way of constructing the new states explicitly is to use the 
following trick. For a given level we solve the null state equation for a state 
with momentum p. Doing so we also obtain a relation between p and the 
central charge c. Rather than choosing a p such that c = 1 and obtaining a 
state which vanishes identically, we keep c ^ 1. Thereby we can isolate the 
zero and extract the new primary. Let us illustrate the procedure by two 
examples. 

The simplest example is obtained by choosing r = s = 1 and p = 0. The 
relevant null state is 

L-i\p >= poc-i\p >, (2.44) 

which vanishes at p = 0. If we divide by p and then set p = we obtain the 
new level 1 primary, 

lim I p >= a_i I > . (2.45) 

A slightly less trivial example is r = 2, s = 1 with p = ^ (a' = 2). The null 
state is 

(L_ 2 - — 1—L 2 _ x )\p >= (p 2 - \)^—{a- 2 - -«_!«_!) |p > (2.46) 
2ytr + 1) 2 p z + 1 p 
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and the new level 2 primaries are given by 



p % ^ (L - 2 -2PTT) L2 l)b >= ± ^(«-^V2«_ 1 «_ 1 )|±^ > 

(2.47) 

The extra primary states can also be understood in terms of SU(2) mul- 
tiplets. If the c = 1 CFT is compactified on a circle with radius R then 
the momentum will be quantized in units of 1/R. When R is equal to self 
dual radius, where the theory is invariant under R —>■ a'/R, there exists an 
extra SU(2) symmetry. The allowed values of the momentum are precisely 
the ones discussed above, namely integer multiples of 1/a/o 7 . The primary 
states, with a given conformal dimension, will arrange themselves into SU(2) 
multiplets. Thus in addition to the discrete momentum tachyon states there 
will exist, for these values of the momentum, additional states which fill these 
out to full SU(2) multiplets. These states are primary then for any circle for 
which the momentum is allowed, and certainly for the real line for which all 
momenta are allowed. A state (r, s), with conformal dimension according to 
( |2.43|) , belongs to a multiplet of dimension r + s + 1. This enables one to 



derive a systematic construction of the new primary states. The states of 
the n + 1 dimensional SU(2) multiplet may be constructed by acting on the 
highest weight state with the SU(2) lowering operator. The highest weight 

state is the tachyon e^ X ^ and the lowering operator is e~ ~7^ x ^ z \ As an 
example, the next to highest state is given by 

dw __2i_xftu~i -^-x(z) f dw : e^ 7 vV : 
• e vV y ' :: e^ 7 



2iri I 2iri (z — w) n 



n-l 



= I 1 ■ e ^ X{z) ~^ X{w) ■ I (2 48) 

Equations ( |2.45| ) and ( |2.47| ), with n = 2 and n = 3 respectively, are easily 
verified, n = 2 gives : dX(z) :, corresponding to the state a_i|0 >. n = 3 
gives 

: d 2 X(z) - 4=(dX(z)f 



'a' 

corresponding to the state (a_ 2 — v / 2a-i«-i)|l/\/a 7 >• 



eVJ X(z) : |0 >, (2.49) 
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When we couple the CFT to quantum gravity the states will be gravi- 
tationally dressed, [pi, p7| 52, 68]. This means that a state with conformal 



dimension h is multiplied by a factor where (f>L is the Liouville field and 
(3 is such that the overall conformal dimension is one. The new gravitational 
scaling dimension and f3 are given by the well known formulae, 




d=l-^ V \ -. (2.50) 

where Q is the background charge, Q = 2\[2 for c = 1. The expression for 
d just comes from comparing the scaling of the specific operator with the 
scaling of the metric according to d = 1 — From this we get the dressed 

dimensions d = f or the tachyon primaries and d = ^ for the null 
descendants, of course still with the same momenta \foHp = r — s. In the 
case of zero momentum which will be our main concern chapter 7, r — s, 
and we find d = r for our new zero momentum primary states. These are 
also the dimensions found in []37| by considering the correlation functions of 
zero momentum operators on the sphere using the matrix model. 

As already mentioned, the special states can be classified using their 
SU (2) quantum numbers J and m. For example, the dilaton is J = 1, m = 
and the state in ( |2.47|) , J = 3/2, m = 1/2. In figure 2.1 a schematic picture 



of all special states is provided. The states at the edges with m = ± J are the 
special tachyons. The SU{2) algebra clearly only connects states with the 
same J quantum numbers. An example was given above in ( |2.48 ). However, 



there exists a much larger algebra , the H^oo, which connects all states and 



of which SU{2) is only a small subgroup. This was shown in [[48] and [|7£ 
In P| with similar methods as above for SU(2). 

As we will see in later chapters there is also a algebra in the matrix 
model. This we will use extensively in some explicit calculations in later 
chapters. 



2.6 Summary 

In this chapter we have looked at some basic properties of the two dimensional 
noncritical string. We have investigated its spectrum discovering the special 
states. We have also considered some simple correlation functions on the 
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Figure 2.1: The special states. 

sphere and in that context noted the peculiar conservation laws for Liouville 
momentum. In the next chapter we will take a natural next step by making 
some calculations not on the sphere but on the torus. 
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Chapter 3 

The Genus 1 Partition Function 



3.1 Introduction 

In this chapter we will use continuum methods to calculate the genus one 
partition function for some low dimensional string theories. Genus one is 
particularly simple since there are no need for any screening charges. We 
will limit ourselves to c = and c = 1. From the conformal field theory 
point of view the first one corresponds to pure gravity while the second one 
describes gravity coupled to a single bosonic field. In string theory we would 
say that we had a one dimensional space time in the first case and a two 
dimensional space time in the second case. 

In || similar calculations were done for any minimal model. We will 
however use a more physical approach and also clarify a few points regarding 
c = 0, "The Theory of Nothing". 

While doing these calculations we will encounter some quite nontrivial 
integrals involving rj functions. In the past such objects have occurred in the 
theory for the critical bosonic string. Due to the presence of the tachyon 
however, the particular integrals have all been divergent and therefore no 
explicit evaluations have been possible. For c < 1 the corresponding integrals 
are convergent and hence it is possible to really do all calculations. We 
will therefore provide some details and examples of such calculations in two 
appendices. 

One motivation for doing a calculation like this is to compare it with the 
matrix model, we will do so and find perfect agreement. Details of the matrix 
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model calculations can be found in a later chapter. 



3.2 String Regulation 

Let us consider a scalar particle theory in d dimensions. Starting with the 
point particle theory we will construct a theory of strings. We will calculate 
the one loop contribution to the vacuum energy. Later we will be interested 
in d = 1 and d = 2 where we may compare with calculations using matrix 
models. This is possible since what we are calculating is, in fact, just the 
genus 1 partition function for two dimensional gravity. 

The one loop contribution to the vacuum energy is given by 

r o V V V r d d v 

E x = - log J D<Pe~ s = - log det □ = -Tr log □ = - J — ^ log □ (3.1) 

where □ = p 2 + m 2 is the inverse propagator and V the volume of space. 

We will now regard this as the low energy limit of a string theory. The 
embedding of the string world sheet in space time we take to be described by 
matter fields of central charge c and also, for noncritical c, the Liouville mode. 
The Liouville mode will appear as an extra dimension, K9, [55]. The volume 



will be the length of space in the Liouville direction. As we have seen, the 
volume is finite due to interactions governed by the cosmological constant and 
given by ( |2.40| ). There are several problems with this expression. One is that 
the momentum integral is obviously UV-divergent. Another problem is that 
it is positive while the matrix models, as we will see in a later chapter, give a 
negative answer. All of this is resolved by recalling that we actually are doing 
string theory. By giving the particle a finite size the integral is effectively 
cut off through the new symmetry which appears: modular invariance. 
If we rewrite: 

V V d 1 r°° 

E x = -Tr log □ = -lim— {— — -/ dtr^Tre^} (3.2) 
2 6 2 s^ods l T(-s) Jo 

t can be thought of as the length around the loop. Since we are dealing 
with strings, the trace (and □) will include higher excited states. For c < 1 
there are no extra field degrees of freedom corresponding to such excitations. 
As we will see below in the case of c = 1, they cancel in the path integral. 
This is consistent with a naive light cone argument where, for d = 2, there 
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indeed are no transverse degrees of freedom. However, as has already been 
explained, there are some special states at discrete momenta. 

As usual we must assure L — L = for the excited states and therefore: 

E , V hm - f 1 [ [ J^ e -t( P 2 +m i) Tr , -tn^(L -L u 

1_ 2 s™ds i T(-s)J V+i2nJ (2n) d 1 



2 s-o 9s T(-s) 7 (vra'ra) 5 ^ 7 (2?r) 

oo 

X e -f (2-rf)-2| J|( 1 _ g n^ 1 _^|2(2-d)| (3 3) 

n=l 

where q = e 2mT . We have renamed <j) — 2hti and t = na'T2 where r = T\+ir2. 
We have used that L + L — 2 is the closed string propagator. The only 
difference from a particle is that we have several excitations to sum over. 

— / 2 

The normalization is fixed by Lq + Lq = + .... The necessary background 
charges, both for matter (if c < 1) and for gravity, only contribute to the 
zero mode parts of L and L , as is clear from ( J2.4jJ . They are responsible 
for giving an m 2 such that modular invariance is obtained. Also, the infinite 
products are the usual contributions from excited states, ghosts giving the 
exponent 2 and matter + gravity the —d. 

The region of integration is originally the infinite strip { — 1/2 < Rer < 
1/2; < Irrir < oo}. The integral is, however, modular invariant. Fixing this 
invariance may be done by choosing just one fundamental region to integrate 
over. The usual choice is F = { — 1/2 < Rer < 1/2; |r| > 1}. This will make 
the integrals finite. We may therefore take the s-limit first: 

V [ £t f d d p ..^y- , , |2CM) 



This term is the one coming from taking the derivative of r(— s), hence the 
minus sign. The other term, from (na'T2) s , is proportional to s and hence 
zero for any finite cutoff. We finally obtain 



23 



3.3 A Theory of Something 



For c = 1 we have d = 2 and therefore according to (|3.5|) 

V f d 2 r V 

E 1 = — / ^ = —. (3.6) 

8ttV Jf r| 247ra' v ; 

Now a' = (a^a^) 1 / 2 , i.e. a' m for matter and for gravity. We also have, 
implicitly, V oc (a^) 1 ^ 2 to get correct dimensions. To be able to compare 
with the matrix model results we put a' m — 1/4, a'* = 2 and use the result 

(EH 

V = --^logA (3.7) 
where A is the cosmological constant to get 

£ 1 = -LlogA (3.8) 

247T 

in agreement with |39|] apart from a factor 1/2 due to doubling of the free 
energy not accounted for there, 0. This will be explained at the end of the 
next chapter. 

This is probably the simplest example of how strings regulate a theory. As 
already stated, the only physical state is a massless tachyon. Furthermore, 
it's only the tachyon which propagates around our loop. Contrary to the 
critical string no excited states are needed to provide modular invariance. 
We have a string regulation of a massless particle in two dimensions. 

This gives a possibility to give a hand waving argument for introducing 
the string. The invariance which turns into modular invariance is inversion 
of the length around the loop t — ► a'/t. We can consistently omit t < (a') 1 / 2 . 
To get this invariance we must, however, introduce an extra parameter <p 
describing the particle. It must transform as — > %-<p to give invariance. The 
consistent way of putting this together is modular invariance. The omission 
of small t is identical to restricting the integration to the fundamental region 
introduced previously. 



3.4 rj Function Integrations 

This section is a mathematical interlude where we will consider the integral 
( |3.5| ) in more detail. The trivial case d = 2, where there are no rj function 
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integrations we have already done above. What about other cases? Below 
we will show how to do the integral for d — 1 and d — — 1. 
Let us first consider d = 1. Our integral is then 

J F $J-2\V(1)\ 2 - (3-9) 
The trick is to rewrite the rj functions in the following way 

\v(d)\ 2 = ^(^E^ |n+H2 - Vf Ee^ |n+mTp ). (3-10) 

V 2 n.m V ^ n.m 



This is proven in appendix 3A. Hence we need to consider 



d 2 T 



£ / ^_ e -^\^\\ (3.11) 



The next step is to use the method of [|5(], |7T| used in ||. The n's and m's in 
( 3. 11|) may be thought of as describing different windings around the torus. 
Since ( |3.11| ) is modular invariant we can use modular transformations to put 
all the winding in each case around one specific cycle. We effectively trade one 
of the sums for a sum over modular transformations. These transformations, 
when operating on the fundamental region, cover the full strip { — 1/2 < 
Rer < 1/2; < Imr < oo}. We get 

d 2 T 



E 



e ^ T 2 

T 2 



2m W+m-rV 



d 2 T f°° dr 2 ^ zh„2 n 2n 

— + 2 / -^E e2T2 =o+tt. 3.12 
f r£ Jo t£ ^ 3 3x 



Putting everything together then gives 

The next example d = — 1 is slightly more involved but the techniques 
are the same. As shown in appendix 3B it is again possible to rewrite the rj 
functions in such a way that the result (|3.12|) can be used. The result is 

d 2 T 



V(Q)\ -TTS- (3-14) 
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The details of the calculation may be found in appendix 3B. 

In this section we have seen some examples of explicitly calculated 77 func- 
tion integrations. Presumably there are other cases where a similar method 
would work. 

3.5 A Theory of Nothing 

Let us now use the result of the d = 1 calculation to obtain the vacuum 
energy for the c = model, pure gravity. The volume in this case equals 
— log A as obtained from Q2.40| ). Note that we have no factor of a' m here. 
The result is 

^ — llogAl, (3.15) 

which again agrees with the matrix models apart from a factor 1/2. 

Let us compare with the derivation in || for c = i.e. the (p, q) = 
(2, 3) model. There a matter field was included and E\ for finite radius R 
calculated. E\{R/ \fc^ = \/pq) ~~ E\(R/ '\/oc' m = \fpjq) then give the (p,q) 
models following f25| , in particular the (2, 3) one. In these cases the Dedekind 
//-functions cancel since we always calculate in d = 2. The complicated 
integrals are instead due to having a finite R. 

For c = we have an alternative, as shown by our calculation above. 
We may do without the matter field completely which is really the natural 
thing for c = 0. In that case the Dedekind functions do not cancel ( |3.9|) , on 
the other hand there are no finite radius to worry about. The equivalence is 
illustrated by the identity ( ft.lOQ . 

So, for c = we may either have no matter at all, really pure gravity (our 
version) or we may have a bosonic field with a background charge giving a 
net c = (as in 0). The two approaches are completely equivalent, and the 
physical content the same. 

In other words, two variants of a theory of nothing. 

3.6 Summary 

In this chapter we have calculated the genus 1 partition function for c = 
and c = 1. We have compared the answers with the matrix model and found 
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agreement. These calculations are basically the only ones which have been 
done at genus 1. This illustrates the difficulties associated with calculations 
using field theory. The more remarkable are the achievements of the matrix 
model to which we turn in the next chapter. 
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Appendix 3A 

We shall first prove 

\ v ( q )\* = W [gf(f+*) 2 gf(f-*) 2 - g f(f+!) 2 gf(f-f) 2 ] . (3.16) 
^ a ,t 

First we note that 

oo oo 

v ( q ) = g V24 JJ(i _ g») = £(-l)" g §(«-*) a . (3.17) 

n=l — oo 

This may be obtained using the Jacobi triple product: 

oo oo 

l[{l-x n ){l + x n - 1 / 2 y){l + x n - 1/2 y- 1 ) = $>* n V, (3.18) 



n=l 



setting x = q 3 and y = —q ^ 2 J3l[. We then get 

^(g)! 2 = ^(_l)n+m g |(n-l/6)^|(m-l/6)^ (3 lg) 

We will prove (|3.16|) directly from this by identifying terms. We first look 

(3.20) 



for terms in (|3.17 ) such that 

(| + t) 2 = (n-l/6) 2 



We find 

{j. n—m 
s = 6{n + m) — 1 

or 

m— n 

s = -3(n + m) + l (3 ' 22) 

with n ± m even. 

So we have found all terms with n ± m even in (|3.17| ) represented. The 
factor 1/2 in (BTTH) corrects for the existence of two solutions above. 
Analogously we find 

( 1 3(n+m) — 1 

\ 2 (3.23) 

s = n — m 
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or 

C , _ -3(n+m)+l 

< 2 (3.24) 

[ s = m — n 

for n ± m odd. 

All terms in fl3.17|) are now accounted for. However, in the first term 
of ( |3.16| ) we have only taken care of terms where s is not even or divisible 
by 3 (see ( |3.21| , |3.22|) ) and in the second term, terms where s is odd and 
t not divisible by 3. What about the rest? They can be seen to cancel by 
considering the solutions to 

\2 _ (82 , t 2 \2 



£2 _ 

2 3 



n _ ^ ( 2 _ ) S2 _ V\2 (3.25) 



i.e. 





= s 2 


t-2 


= 3ti 




= 2t 2 


s 2 


= 2ti 



(3.26) 



or 

f c — 9f„ 

(3.27) 



(|3.16| ) is now proven. Finally the formula for Poisson resummation is needed 



1 



CO 



n=— oo m=— oo 

see e.g. f35j. A resummation on s yields the desired result (|3.10| ). 
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Appendix 3B 

To calculate 



' r ' T 'vW)f (3.29) 



F 1/2 
r 2 



we again appeal to the Jacobi triple identity as in appendix 3A. In this case 
it yields 

}J(l-q n f= hm ^-°° g i/2 ; . (3.30) 

1 /2 

With y = —jz^ this becomes 



1 / -1 \ n 



n(i-g n ) 3 = iimE-(T^) <P n ^ n 

n=l £ -00 e \ 1 e / 

and finally 

00 

^(g) 3 = ^(2n + 1)(-I) n g5(«+D 2 . (3.32) 

n=l 

We then need to prove 

I 77(g) | 6 =4 J2 {n+-){m+-){-l) n+m q^ n+1 2) 2 q^+\? 



(3.31) 



n,m=l 

00 2 

= 2 E (j-^qKWqto-*?. (3.33) 

s=odd,t 

This can be proven in a similar way as in appendix 3A by identifying terms. 
Then we put q = e' l7rT and obtain 

00 2 2 
I v (q) | 6 = 2 2 ( j - ^e-^tV+^w 



00 , 2 1 J 00 

2E<=-'"-(T + if-)Ee 



27rr 2 t 2 +2i7rrist 

4 ' 2n dr 2 '' 



s=odd ' 1 t 



2 y g-WTa-y _ + )[^ =e 2r 2 y e 2r 2 r 2 3.34 

s^dd l 4 2WW^ 4- 
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where we in the last step performed a Poisson resummation. 
From this we find 



F T 1/2 



1 <~-b E /$^ IW " P U, (3.35) 



7i V2 dx 2 s J^ J t. 
Let us calculate 



£ J d -L e -^\^ (3.36) 

s=odd,n 2 

The only slight complication is that the sum is only over odd s. Under 
modular transformations the measure is invariant and in the exponent we 
see that 

t — > t + 1 takes s — > s,n — > n + s 

t — > — 1/t takes s — > n, n — > — s. (3.37) 

We therefore realize 

E = E = E • (3-38) 

s=odd,n=even s,n=odd s=even,n=odd 

We can not, however, get J2s,n=even i n this way. But this sum is immediately 
given by ( |3.12D to be 

e (3-) 



s,n=even 



We therefore find 



c/ 2 r 
r.- 



7— In+srH 
g 2t 2 1 1 



E + E =~(E- E ) = T7T- (3-40) 

s=odd,?i=et)en s,n=odd ° s,n s,n=euen OJ ' 



By insertion in ( |3.35| ) we then get 

d 2 T 
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Chapter 4 

The Matrix Model 



4.1 Introduction 

Matrix models are new and powerful tools which can be used to solve a 
variety of different models for two dimensional quantum gravity. Hence they 
provide solutions of noncritical string theory. To be honest, matrix models 
have been around for quite some time as models for random surfaces and 
two dimensional quantum gravity. In the beginning however, one usually 
limited oneself to a large N limit (N being the dimension of the matrix) such 
that only surfaces of spherical topology survived. It was a very important 
discovery when it was first realized that an appropriate continuum limit, "the 
double scaling limit" , existed which made it possible to retain all genus in a 



sensible fashion. The first such matrix models || £3 |1[ described various 
theories of quantum gravity coupled to, not necessarily unitary, matter with 
central charge c < 1. The model which will be of interest to us, however, is 
the case c = 1, the topic of this thesis. While c < 1 was basically a problem 
of doing very complicated matrix integrals, something which may be done 
using orthogonal polynomials, the c = 1 case requires a different method 
In this chapter we will review some of these issues. 
Unfortunately the applicability of the matrix model seems to stop at 
c = 1. This is similar to the techniques of field theory in previous chapters. 
It is presumably related to the fact that for higher values of c one gets new 
field degrees of freedom which cannot be encoded in just a matrix. We should 
also recall that the presence of a tachyon makes the d > 2 theories quite ill 
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behaved. A natural way to proceed would therefore be to consider super 
symmetric theories. When this is being written, however, the progress in 
such directions has been very limited. 

4.2 The Matrix Model 

Let us start with the c = 1 string partition function 

Z = ± 9T 1] I VhVXe^ ^ (4.1) 

9=0 J 

where g s t is the string coupling. We have summed over all genus g. In the 
path integral we have included an integration over h, the two dimensional 
world sheet metric. Classically, the action is of course independent of h due 
to Weyl and reparametrization invariance. Quantum mechanically, as we 
have seen, the measure VX has a dependence due to the need of a regulator. 

The basic idea of the matrix model approach is to represent the surface 
by a triangulation. The triangulation, or rather its dual, is regarded as a 
Feynman diagram of some scalar field theory. The reason to consider the 
dual graph is of course that the order of all the vertices is then the same. 
It is three in the case of a triangulation although any polygon would do. 
The vertices correspond to interactions and the edges to propagators in the 
scalar field theory. The matter, i.e. X field, part of the action becomes in 
this discretized version 

j VXe -^^ w ~/n^^i' (4.2) 

i 

rf is ~ I/a', < ij > means a sum over nearest neighbors. 

In this picture the integration over all metrics, / T>h, turns into a sum 
over all triangulations. From ( J4.2p it is clear that we need a space time theory 
whose propagator is 

e ~v 2 \u-t 3 \\ ( 4 _ 3 ) 

1 p 2 

Since the Fourier transform is given by e *r , this requires the very com- 

_ 1 d 2 

plicated kinetical term <pe Such a theory would be very difficult to 

handle. Fortunately we do not need this. In the continuum limit the lattice 
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spacing is going to zero, and the space time momentum is going to zero in 
units of one over the lattice spacing. In the continuum limit we are indeed 
not supposed to see the discreteness of the surface. This is precisely what 

1 -2 

we need. In this limit the propagators e **r and p2 ^ 4r)2 coincide. It may 
be noted that , }. 2 is the Fourier transform of e _T? l* i_tj 'L Hence the model 
to investigate in the continuum limit is simply 

2 + V0 2 + V((j>) (4.4) 

where we have added some interactions V(4>). A 3 would give triangula- 
tions. The argument above should not be thought of as a proof. It is rather 
a heuristic argument which serves as a motivation for further study. The 
models ultimate success and agreement with known results from Liouville 
theory is the best justification. 

Given a triangulation it is a classical result that the genus g of the corre- 
sponding surface is given by 

2-2g = V-E + F (4.5) 

where V is the number of vertices, E the number of edges and F the number 
of faces of the triangulation. It is important that our scalar field theory can 
measure the genus, otherwise we could not introduce the string coupling g st . 
The correct way to achieve this is by letting <fi be not just a single scalar field, 
but an N x N hermitean matrix. The partition function to study is then 

Z = J V(f>e-^ Tr ^ 2+u ^\ (4.6) 

where (3 is 1/h. From here we see that each face of the triangulation is 
weighted by j3, from the Feynman diagram vertex, each edge by (3 from the 
propagator, and finally each vertex or Feynman diagram loop by N. This 
last point is the reason for having matrices. Combining this we find 

{±) E P F N V = (3 2 - 2 °(j) v (4.7) 

for the weight factor of some given triangulation. Clearly the string coupling 
will be g st ~ 1/ [3. If this is to describe a continuum string theory there must 
be a continuum limit such that surfaces with infinite number of vertices, or 
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faces, dominate. From above we understand that the partition function can 
be seen as a power series in ^. The trick is to adjust ^ to some critical value 
where this sum diverges. This means that the V — > oo surfaces dominate 
and we have a continuum limit. However, we clearly want to extract some 
sensible, finite answers. That this is possible is the single most remarkable 
feature of the matrix model. By taking iV and (3 to infinity in such a way that 
the critical point is approached, the resulting divergence may be cancelled 
against the 1/(3 behavior of the string coupling, renormalizing everything to 
finite values. This is the double scaling limit. It is a nontrivial fact that this 
is possible to do simultaneously for all genus. 

Now that we have established this equivalence of models, how do we solve 
the matrix model? A hermitean matrix <p may be decomposed into a diagonal 
piece with eigenvalues and a unitary part U, i.e. = UAW. The matrix 
model measure is similarly decomposed into these parts 

V(f) = Y[dX i VUA 2 (X i ) (4.8) 

Xi are the eigenvalues, and 

A(A0 = II(Ai - A,-) (4.9) 

i<j 

the Jacobian of the change of variables. It is sometimes called the van der 
Mond determinant. 

If we integrate over the unitary part U the remaining partition function 

is 

Z = j PAA(A i (t 2 ))e" /3 ^i 2 i^?+tf(*0] A(A i (t 1 )), (4.10) 

where the eigenvalues Xi are fermionic. This procedure is strictly speaking 
only valid for the infinite real line (for the t variable). At finite radius the 
angular variables do not decouple. We will come back to the finite radius 
case in a later chapter. 

The double scaling procedure is controlled by introducing the critical 
parameter 

A-l-£, (4.11) 

assuming that we have renormalized in such a way that ^ = 1 is the critical 
point. Since V measures the area of the world sheet we understand from 

{!L ) v = ( i- A) v„ e -*y (4.i2) 
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Figure 4.1: The top of the matrix model potential and the Fermi sea. 



that A must be the world sheet cosmo logical constant. 

It is now time to attempt some more explicit calculations. A suitable first 
object is the ground state energy E = E(A) of the system. First we note 
that a small change in A is the same as a change in the number of fermions 
N. Clearly, in the ground state, the fermions are filling up the potential well 
up to some Fermi energy \x. Therefore 

^ = fr (4.13) 

We have chosen the zero of the energy to be at the top of the potential, see 
figure 4.1. 

Since n has a very clear matrix model interpretation, it is more natural to 
work with \x than with A. In this picture the double scaling limit is obtained 
by taking (3 —>■ oo while keeping nf3 fixed. In this limit the only thing that 
matters is the top of the potential. The top is magnified more and more as 
(3 — ► oo and fi —>■ 0. Even so, it is clear that some kind of a cutoff is needed 
for the proper definition of the theory. The Fermi sea must have a bottom 
and another shore somewhere. The typical case is a potential of the form 

H 

^<A'-A<). (4.14) 
When we expand around the top we find the leading contribution 

We have redefined A by a shift to the maximum of the potential. Another 
type of cutoff which also is relevant for a precise nonperturbative definition 
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of the theory, is to put infinite walls on both sides of the potential. In either 
case, as the double scaling limit is approached all higher corrections to Q4.15Q 
become irrelevant and the walls move off to infinity. This is important for a 
well defined theory. 

Let us now introduce the Legendre transform T(p) given by 

E(A) = P 2 iiA-T(ii). (4.16) 

By standard arguments E(A) is the generating functional for connected di- 
agrams, while generates 1PI (one particle irreducible) diagrams with 
respect to the puncture. We will, unless stated otherwise, consider the 1PI 
amplitudes. The notation will be 

d n T 

< P...P >= —— (4.17) 

dp n v ' 

It is very convenient at this point to introduce the density of states, p(p), 
clearly given by 

/ s 9A , 
P(H) = ^ (4.18) 

Using this, the ground state energy can be written as 

E(A) = ep(e)de. (4.19) 



Note that this is consistent with 

dE dE dA 



Mm) (4.20) 



dp dA dp 
if we think of A as A(p). Also 

A = p(e)de =< P > . (4.21) 

This is a good point to introduce the Fermi liquid picture which uses phase 
space to illustrate the system. In phase space each fermion follows a trajec- 
tory p 2 — A 2 = e, where e is its energy. Hence the system can be described, 
classically, by figure 4.2. A liquid rotating in time. The excitations, i.e the 
tachyons, correspond to ripples on the surface which rotate, and evolve, with 
time. 
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Figure 4.2: Phase space picture of the Fermi sea. 

With phase space coordinates the puncture one point function can be 
written as 

< P >= [ dpd\. (4.22) 
The two point function follows from this as 

< PP > =—J dpd\Q{-n - H) 
= I dpd\5(-i2 -H)=l . (4.23) 

J J Fermisurface 

More generally for some operator O consisting of polynomials in A and p: 



<OP>= f O. (4.24) 

J Fermisurface 



The two point, or p(fj), can also be written as 
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< PP >= -TmJ2 



it n e n + t + it 

I /*oo j fOC f 

= -Im / dTV < n | e~ T ^ +H) \ n >= -Im / dT d\G(X,X;T). 
n Jo n n Jo J 

(4.25) 

We see that there are several ways of approaching computations of matrix 
model correlation functions. One way is to start with the first expression 
in (|4.25|) . One may then analytically continue to a right side up oscillator. 
This was the way the original calculation, [42]j, was done. The sum over the 



energy eigenvalues becomes a sum over the imaginary energy eigenvalues in a 
continued right side up. This is also the way in which we will do computations 
in this thesis. Another approach is to use the path integral formulation, i.e. 
to use the last representation in ( |4.25| ). This was, for nonzero momentum, 



first done in |Sq] and extended in 

For reference we give the genus expansion of the puncture two point 
function as obtained in [E2] 



< PP >= --Im Y = -l^RetP C ^ v "" u ) 

~ ^ 0i ^(2n + l)-t tt n 2 



7T 

where B 2n {\) = {2 l ~ 2n - l)B 2n . 

This is a good point to reconsider the results of the previous chapter. 
From the relations derived earlier in this chapter, it follows that for genus 
one 

S 1 = -_J_|logA|. (4.27) 



24ttvV 

As we have pointed out this is twice as large as the field theory answer. The 
explanation is however simple. Recall the matrix model potential and the 
Fermi sea. The expression ( f4.26|) really presupposes that there are two Fermi 
seas. One on each side of the top of the potential. Hence there are really 
two copies of the world, perturbatively disconnected. Clearly this leads to 
a doubling of the free energy. In the next chapter we will briefly consider a 
case where there is only one world. 
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There are several additional comments to be made about the genus 1 
partition functions. This will also, in fact throw some light on how special 
the choice of an inverted harmonic oscillator is. 

As discussed in chapter 3 the genus 1 partition function is basically ob- 
tained just from the diagram of a single particle closed loop. This we could 
write as ~ / d 2 G(x,x), where G(x,y) = log | x — y | is the propagator at 
equal times, x is a space coordinate. This certainly needs some regulariza- 
tion. One possibility is a lattice cutoff L in some finite volume L. The result 
would schematically be 

Ei ~ L(i + ^) (4-28) 

The second term is the ( function regularized sum of oscillators in the volume 
L. In the L — > oo limit this term is dropped. It is only this second term 
which is universal, the first one needs some new physics, e.g. the scale of 
the cutoff, to be determined. The modular invariance of the last chapter is 
precisely such a prescription. The expression (|4.28|) looks rather symmetric. 
In fact, the first term can also be thought of as a ( function regularized sum, 
but now in some "internal" space of volume L. We will see in a moment from 
where this comes in the matrix model. 



As noted in ||22|| , the E\ can be thought of as arising from insisting that the 
coincident propagator of the one loop diagram should be defined by normal 
ordering in A space. When we coordinate transform to r space, where r is 
the time of flight for A, we pick up a Schwarzian derivative which precisely 
gives Ei. The map between A and r is A = y^coshr in general. At the top 
it is A ~ e T . But this reminds us about a similar situation in conformal field 
theory where one maps from the plane to the cylinder. There it is well known 
that the Schwarzian derivative of the map and an explicit sum over modes 
on the cylinder give the same result. In our case the circumference of the 
"cylinder" is the imaginary period of e T , i.e. the time period for oscillations 
in the continued harmonic oscillator. This is L. As we will see later, this 
period is also related to the quantized momentum of the special states. 

This discussion began with a promise that it would say something about 
the more general models with anharmonic potentials. It is clear that L, the 
characteristic size of the string, is independent of the cosmological constant, 
either fi or A, only for the harmonic potential. Hence we have L ~ log fi and 
L ~ 1. In general the period depends on the amplitude and we have both L 
and L ~ [i~ k for some positive k. 
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Not only does the effective size of the string depend on the cosmological 
constant, there is no natural way of keeping the scale of the string and the 
scale of space very different. This is only one of many strange features of the 
multi critical models. 



4.3 The Special States 



In the discussion of the Liouville theory in chapter 2, the special states were 
introduced. 

In the matrix model representation of the theory it is easy to see some 
of these degrees of freedom. In particular, there exist special operators of 
zero momentum, which are represented as time independent perturbations 
of the matrix model potential. These are the analogs of the scaling operators 
in the one dimensional matrix model. They have the physical effect, in the 
matrix model, of moving one from one critical point to another. Their cor- 
relation functions are easily computed on the sphere and for a general genus 
surface. In chapter 7 we will investigate these zero momentum correlation 
functions and show that they obey certain recursion relations which relate 
all of them to the puncture two point function. We will consider both one 
particle irreducible and connected amplitudes. The recursion relations allow 
for computation of any zero momentum correlation function for any genus. 
They are very similar to the recursion relations for c < 1 theories and suggest 
a topological interpretation. We will also rewrite the recursion relations in 
the form of constraints on the puncture one point function. These constraints 
obey a Virasoro algebra. 

We will then proceed, in chapter 8, to consider the operators of nonzero 
momentum. The obvious candidate in the matrix model for these operators 
are the time dependent perturbations of the matrix model potential. We will 
calculate correlation functions both on the sphere and on a surface of arbi- 
trary genus of such time dependent perturbations. We find that correlation 
functions of these operators have real poles in the Euclidean momentum, 
conjugate to the target space variable, at precisely the expected values of the 
momenta of the special states. 

What is the meaning of these special operators in the matrix model? 
A hint to their meaning is provided by the analysis, f37|, of the discrete 



string. This is the model in which the string is mapped onto a discrete set 
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of points with equal spacing e. It was shown that this model was equivalent 
to the model where the string is mapped onto a continuous line, as long as 
e < e cr , at which point a transition appears to take place. Now this model 
of triangulated surfaces is the dual of the case where the string is mapped 
onto a circle of radius R = - and the transition is dual to the Kosterlitz- 
Thouless transition which arises due to the condensation of vortices. We will 
discuss this a little bit more in chapter 6. One explanation of the equivalence 
of the continuous and discrete strings, given in fl38| , is to note that the 



discretized line can be replaced by a continuous variable t with a periodic 
potential, V(t) = Xm a n cos(— t). The lowest dimension operator in this 
potential, n — 1, has undressed conformal dimension ^(^r) 2 and becomes 
relevant (i.e. of undressed conformal dimension one), when e = e cr = tt^/cV. 
This is where the phase transition of the discrete string took place. The 
variable t is forced to take values at the minima of V(t) and the real line is 
effectively discretized. Now, these terms in the potential correspond precisely 
to periodic time dependent perturbations of the c = 1 matrix model and the 
place where the transition takes place corresponds precisely to one of the 
special values of the momenta at which the special states occur. Thus we 
might conjecture that the physical meaning of these operators is that they 
transform the continuous line into a set of discrete points, with a period 
given by the inverse of the special values of the momenta — spontaneous 
punctuation. 

As we will see later on, it is not enough to just consider the matrix 
eigenvalues when identifying the special states. To really be able to see the 
individual special states one needs to invoke the conjugate momentum. The 
precise way in which this leads to the symmetry will be discussed in 
chapter 9. 

Finally it is useful to point out that one important advantage of the 



analytical continuation approach advocated in |15|, [T(|, is that the special 
states are given a very clear interpretation. As we will see in more detail 
in chapter 8, they correspond to resonances in the harmonic oscillator. This 
is also a reflection of the fact that they are present for Euclidean momenta. 
Hence to really see them in a Fermi liquid picture, we must use the analytical 
continuation. This is illustrated in figure 4.3. Now any special state is 
associated to some periodic ripple on the Fermi surface. 
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Figure 4.3: Analytical continuation of the phase space picture. 

4.4 Summary 

In this chapter we have described the c = 1 matrix model and how it solves 
many of the problems of formulating a noncritical string theory. Several 
formulas for calculating correlation functions have been obtained which will 
be useful in later chapters. We have also encountered the special states which 
will be of major importance in coming chapters. 
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Chapter 5 

Nonperturbative Effects 



5.1 Introduction 

In this chapter we will discuss some of the nonperturbative issues of two 
dimensional string theory. This is really the most remarkable property of the 
matrix model solution, the possibility of obtaining nonperturbative results. 
A general discussion of nonperturbative effects in string theory |75| suggests 
that the (2g)\ behavior in ( |4.26| ) and the e~ x l 9at behavior we will find below 
are generic for closed string theory. This was first observed for the c = 1 



model in Wl 



5.2 One World or Two 

We will consider an explicit and very simple example as an illustration of 
what happens nonperturbatively. An extensive treatment can be found in 
For our limited purposes we will use other and simpler methods. 
Our starting point, as for most cases in this thesis, will be the expres- 
sion ( |4.2C ) for the puncture two point function. Perturbatively each of the 



terms in the genus expansion can be obtained using the WKB approximation. 
Since the expansion turns out not to be Borel summable, the summation is 
not perturbatively determined. Additional information is needed. Roughly 
speaking, we have the freedom to add terms of the form e~ 1//9s ', where g s t is 
the string coupling. Such terms do not show up in the small g st expansion. 
The two point function as given by ( f4.26| ) therefore includes additional non- 



44 



perturbative information which seems to have been provided by the matrix 
model. However, as we will see, the choice is not unique. 

The following is a simple modification of the matrix model which leads to 
a nonperturbatively different result. Instead of having just the upside down 
harmonic oscillator potential, we imagine having an infinite wall erected at 
the top of the potential. One might say that one has one rather than two 
worlds. Clearly the WKB expansion, which only feels the potential up to the 
classical turning point, will be insensitive to this. The wall can only be felt 
nonperturbatively. The expression for the puncture two point function can be 
easily adjusted to describe this situation. The effect of the wall is to impose 
a new boundary condition by forcing the wave functions to be zero at the 
position of the wall. Hence, rather than summing over all energy eigenstates, 
we sum only over the odd ones. Therefore, the two point function in the case 
of the wall is given by 

< PP >= iRe V t - 1 , - = -^Re^( 3 + 2v/ ^ ). 

(5.1) 

By the above reasoning we know that this must perturbatively be precisely 
half of the result without a wall. This can be checked explicitly by using 
the same kind of expansion as in Q4.26|) and the Bernoulli polynomial iden- 
tity £>2n(l/2) = 4 n i?2n(3/4). However, nonperturbatively there may be a 
difference. Let us calculate it. By using the ip function identities: 



■0(1 — z) = 1p{z) + 7T COt TTZ, (5-2) 



one easily shows 



^(—) =^(— )+log2--cot(7r— )+-(^(— ))■ (5-3) 
With z = 2\fa! (3 ixi imaginary it follows 

Re ^( ^^ )=Re ^ ( ^^) + log2- . (5.4) 

^ 4 2cosh7rvapA i 
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The nonperturbative nature of the difference is obvious. Expanding it for 
small e~ l l 9st we find 

oo 

_ n Y, {-l) n e - {2n+l)nV ^^. (5.5) 

n=0 

This nonperturbative ambiguity can also be understood from an integral 
representation point of view. We note that 

^(^) = n 1 -^) + - ^(^)). (5.6) 

Using 

/■oo fp~ z ^ 

^>{z) = / - (5.7) 

Jo 1 — e * 

one can, by a change of variables t — > —t, prove that for z = 2\fa'l3jjii 



-Re^(^) = Re jf 



a' 4 io 1-e" 



l + z , 



R^{—-) = Re / — — . (5.8) 



r a' dfi\i 4 Jo 1 — e 



The integrands are the same, but the contours differ. One is the real line from 
to +oo, the other one the real line from to — oo. If we rotate one contour 
into the other, we pick up residues from the poles along the imaginary axis, 
this precisely gives (|5.5| ). 

Hence we have established the nonperturbative nature of the difference 
between one and two worlds. 



5.3 Instant ons 

It would be nice to have a more physical picture of the nonperturbative string 
theory. We will try to achieve this below. 

Nonperturbative effects are in general associated with instantons. In this 
case we have eigenvalue fermions at the Fermi surface tunneling through the 
upside down potential. Without the wall, the instanton action for one fermion 
tunneling through is obtained from the Euclidean action (3 J <it(|A 2 — ^A 2 ) 
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evaluated for the classical tunneling solution of energy \i. We also need to 
subtract the action for a fermion not doing any tunneling, i.e. just sitting at 
ji — 2^7 A 2 . The result is 



13 / , dX\2/i - —A 2 = Vc/nfip. (5.9) 

V a' 

Now, for two worlds the instantons correspond to tunnelings through 
the potential and back again. The leading instanton contribution to the 
propagator, G(X, A), has then the weight 

e- 2 ^ 7 ^. (5.10) 

Considering the integral representation of Reipi^rr) we indeed note poles 
leading to such contributions. 

Let us turn to the case of one world. We write the two world propagator 
as a sum over even and odd (i.e. one world) parts G(A, A) = G e (A, A) + 
G (A,A). Then G(X, -A) = G e (X, A) - G (A,A) and from that G(X, A) - 
2G (A, A) = G(X, — A). The leading contribution to G(A, —A) is a tunneling 
solution through the potential barrier, but not back again and hence of order 

This verifies and illustrates the calculations in the previous section. 



5.4 Summary 

Clearly much remains to be done concerning the nonperturbative aspects of 
the two dimensional string. The important contribution of the matrix models 
is a laboratory where we can check our understanding of nonperturbative 
string physics. As we have seen above, the nonperturbative effects have a 
rather simple interpretation in terms of single eigenvalue tunneling. What 
remains to be obtained is a clear picture of what happens in a string language. 
One possibility, which has been raised in the past, is that by clarifying this 
issue, one could eventually find the "real meaning" of string theory. Perhaps 
this will involve concepts which transcend the traditional genus expansion. 
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Chapter 6 

String Theory at Finite Radius 



6.1 Introduction 

In this chapter we will consider the two dimensional string theory with com- 
pactified time on some radius R. In other words, string theory at finite 



temperature. This has been extensively treated in [pi]] . We will however use 
another and very simple approach to reproduce some of these results. It is 
our feeling that the method below may provide some useful insights into the 
basics of the matrix model description of string theory. 

As mentioned in chapter 4, it is not really correct to ignore the angular 
parts of the matrix integral for finite R, |37|, p8| . The angular parts are related 
to vortices. These vortices are found to cause a phase transition at a radius 
of twice the self dual radius, this is the Kosterlitz-Thouless transition. As 
described in [[ft], the transition is caused by operators cos(-^(X R — X L )), 
which create vortices and become relevant at R = 1\fo! '. Below this radius 
the vortices destroy the string world sheet. Therefore, the calculation which 
we will consider below, assumes automatically that vortices are excluded for 
some reason. It is not clear if this is a physically reasonable assumption. 

In the presence of vortices the target space duality of string theory under 



R — > a' I R is explicitly broken. One might comment, [64 , that there is a 
dual object to vortices, "spikes", which tend to condense for large R, more 
precisely R > \\fa' . The spikes are created by cos(^-(X/j + Xl)). This is 
the punctuation phenomenon we discussed at the end of chapter 4. If we 
include such spikes, duality is restored but to the price of having an instable 
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theory for any radius! Luckily, spikes may be forbidden simply by imposing 
translational invariance since the operators which excite them, break this 
invariance. This is not true for the vortices. 

The fact remains, however, that duality is broken. Clearly, this could 
have far reaching consequences also for higher dimensional string theories. 
The argument leading to the Kosterlitz-Thouless transition is, after all, a 
very general one. It is just, as usual, a matter of competition between energy 
and entropy. To see this, one has to introduce a cutoff, like the triangulation 
of the matrix model. As the cutoff is removed, the energy of a vortex goes to 
infinity but at the same time the number of possible vortices also rises. The 
temperature decides which one will win. 

Hence we might learn from the matrix model two dimensional string the- 
ory to be careful when we apply concepts like duality to draw general string 
theory conclusions. We might be using a principle not realized in a realistic 
theory. 

With this in mind, let us now give a very simple finite radius calculation. 



6.2 A Simple Example 

Let us derive a finite radius expression for the puncture two point function. 
To do that, let us consider the fermionic field theory action on finite radius 

S = ['\t d ^ + ^f x+ U M (6.1) 

The path integral, 

Z = J V^V^e"^, (6.2) 

is easily done by expanding ip in eigenf unctions of in time and eigen- 
functions of the harmonic oscillator in the A coordinate. The eigenvalues 
are 2w ^ i 1 for the time part since we need anti periodic boundary conditions 
for the fermions. For the harmonic oscillator part we have the standard 
(2n + 1). There is also the energy shift t = P^- A standard evaluation 
of the logarithm of the path integral, giving a logarithm of a determinant 
and hence a trace of a logarithm of the eigenvalues, then yields the following 
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expression for the two point function 

< pp >= T 1 Fo -tv • (6.3) 



If the sum over m is performed the result of |j7[] is recovered. Using this ex- 
pression a general correlation function can be obtained through perturbation 
theory in the same way as for the R = oo case extensively treated elsewhere 
in this thesis. 

Once this is established, it is also easy to write down a recipe which 
converts the ordinary R = oo two point function into the finite R case. The 
point of doing so is that this procedure, as is easily seen, commutes with the 
variations giving other correlation functions. As we will see, this procedure 



is precisely the one derived independently in |5T 

i, what one has to do is to take one u 

~~2R 

a 



In words, what one has to do is to take one to derivative and then sum 
r all shift 
the identity 



over all shifts ^t 1 ^ j n the to variable. The shift can be written down using 



e z ^f(x) = f(x + z) (6.4) 
for any function f(x). We get 

id 1 1 

(1— e Rat °)y- : — — — — (2m+ i)i 777 = y, 



n,m 



■>FI 



(^(2n + 1). + ^ - t ) 2 V (^(2n + l)i + ^ - t ) 

§- T-^ 1 __L^_ d t-^ 1 

'*o \ — , = e 2R 9t o — > — i 

» ^^ 2n + ^ ~ *°) 2 dto n^r(2n + l)i - t 

(6.5) 



Hence 



i a 



J_ \p 1 _ Rdtg V - 

R fa (^(2n + l)i + - tQ )2 ~ _ e -^4 % ^=(2n + l)i - t ' 

(6.6) 

This verifies the result of BT[| using a different method. 



There are some interesting features of this derivation. Duality is very 
explicit in the expression ( |6.3D . That is, if we do not care about the vortices 
as discussed in the introduction. We can in fact learn a bit about how it 
is realized in the matrix model. Apparently duality relies on the fact that 
the eigenvalue spectrum of a harmonic oscillator is the same as for a fermion 
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with anti periodic boundary conditions in a box! Clearly one cannot expect 
duality to hold in any matrix model. In general a nonharmonic potential 
would have nonequally spaced eigenvalues and hence break duality. 

Curiously, by the same reasoning, many correlation functions will break 
duality. This is not surprising for nonzero momentum correlation functions. 
However, the above expressions suggest this to be the case even at zero 
momentum for the special states. 

6.3 Summary 

In this chapter a simplified derivation of some results for finite radius has 
been given. We have indicated how duality is realized in the matrix model. 
An important observation is that naive duality is not expected to hold in a 
general matrix model. 
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Chapter 7 



Zero Momentum Correlation 
Functions 

7.1 Introduction 



A first exercise in the two dimensional string theory is the calculation of 
zero momentum correlation functions. Clearly this is the simplest thing to 
do. Such calculations will be described in this chapter. We will, apart from 
some comments towards the end of the chapter, limit ourselves to powers of 
the matrix model eigenvalue. The conjugate momentum will be introduced 
in later chapters. One of the main results of this chapter will be a set of 
recursion relations for these correlation functions. 



Let us briefly recall the basics of the matrix model approach to c = 1 quantum 
gravity or two dimensional string theory as described in chapter 4. We will 
limit ourselves to the uncompactified case. The partition function is given 



where is a hermitean N x N matrix. As shown in , the angular variables 
may be integrated out to leave only the N eigenvalues. They correspond to 



7.2 The Gelfand-Dikii Equation 



by 




(7.1) 
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N fermions in the potential U. Below we will consider correlation functions 
of precisely these eigenvalues. 

As we have seen, the density of states, p, in the potential, is the same as 
the 1PI puncture two point function and hence a very interesting and basic 
quantity in the theory. In the coordinate basis it can be written as 

If 1 

p(p F ) = —1m / dX < A | | A > (7.2) 

Trp J H - p F - ie 

By investigating some properties of this expression, we will be able to ex- 
tract information about zero momentum correlation functions of the matrix 
eigenvalue variable. 

We begin by defining the function 

*W =< " l -i v ;Ua) 1 A >' (7 ' 3) 

i.e. the diagonal component of the resolvent. This function obeys the non- 
linear Gelfand-Dikii equation 



- 2RR" + {R'y + 8U{X)R 2 = 4. (7.4) 

Here A is shifted so that the maximum of the potential is at the origin. We 
have also shifted U by pp = p c — p. In the double scaling limit we expand 
U(X) around its maximum and rescale the coordinate. If U(X) ~ A n , then 

2 

A ~ (3 n + 2 . In the following it is assumed that all quantities are rescaled in 
the appropriate fashion. 

The correlation functions we are going to calculate will be correlation 
functions of time independent perturbations, A fc , around the top of the po- 
tential. The corresponding operators we call Ok- They are essentially the 
operators ^Tr[0 — (p c ] k . To perform the calculations we include in U general 
perturbations X k coupled to sources tk, which we later will vary. Taking one 



derivative of (7.4) gives the linearized version of the Gelfand-Dikii equation, 



- -R'" + AUR' + 2U'R = 0. (7.5) 

This formula is also easily verified by using R(X) = J2^=o ^ n ^ x ^ n ^ , and the 
Schrodinger equation for the energy eigenfunctions tp n (X). Note that we are 
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dealing with an inverted potential and the E^s are therefore, in general, 
complex. It is this linearized version of the Gelfand-Dikii equation which we 
will exploit in the following in order to get information about the correlation 
functions of the theory. 

The resolvent R(\) determines the density of states, 

p(p) = -Im / R(X)d\, (7.6) 

7T J 

which is also the puncture two point function < PP >. Other two point 
functions can be obtained by taking the moments of R, 

< O k P >= -Im / X k R(X)dX. (7.7) 

7T J 

These amplitudes are one particle irreducible (1PI) with respect to the punc- 
ture. 

Some comments on the integration region are in order. There are essen- 
tially two different cases depending on whether we have one or two worlds as 
explained in chapter 4. From the point of view of the WKB approximation 
we should integrate only up to the top of of the potential, i.e. from — oo 
to zero in the case of one world. With two worlds, we should integrate all 
the way from — oo to +oo. This is effectively a doubling of the system. At 
a critical point, where we just have a leading even power contribution to 
U(X), R(X) is clearly even. This immediately implies that, with two worlds, 
( |7. Tp is zero for k odd. For the quadratic critical point this is true also with 
one world. The above integral for an odd operator turns out to be real and 
therefore gives no contribution to the imaginary part. 



7.3 1PI Recursion Relations 

We will now use the linearized Gelfand-Dikii equation to derive recursion 
relations for the 1PI amplitudes. We start with ([7.5|) for a general potential 



U(X) = Y,t P X p . (7.8) 

p=0 

The tp's will act as sources for the O p operators. We will later restrict 
ourselves to the quadratic critical point 

U = V^0,2 t = pfi t 2 = __L<o. (7.9) 
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We then compute the A fc+1 moment of ( [7.5|) , obtaining 



J oo 

-~(A; + 1)A;(A;-1) < O fe _ 2 P > +2^t p (2A; + p + 2) < O fe+p P >= 0. (7.10) 

To obtain this, we have made repeated partial integrations and discarded 
boundary terms at infinity which do not contribute to the singular part of 
the correlation functions. In the case of one rather than two worlds, one 
might worry about boundary terms at A = 0. However these are real, as 



follows immediately from the WKB expansion of R(X) in terms of yU(X). 
Hence, they do not contribute to the imaginary parts of the integrals, which 
give the correlation functions. Using < O p ... >= ■£-<••■> we may take t p 
(p > 0) derivatives of ( |7.10| ) to derive 

-\{k - l)k(k + 1) < O k - 2 O kl . . . O kp _ L P > 
+4t (k + 1) < O k O kl . . . O kp _,P > +U 2 {k + 2) < O k+2 O kl . . . O kv _,P > 
+2E|Ti(2A; + h + 2) < O k+k p kl ...6 kv ..O kp _ 1 P >= 0. 

(7.11) 

For simplicity we have chosen the quadratic critical point of (|7.9| ). When 
we take the derivatives, t = (3^ is kept constant. This is because (5\x is 
the nonamputated puncture one point function and we are considering 1PI 
amplitudes. 

With this recursion relation any zero momentum correlation function for 
arbitrary genus may be expressed in terms of the puncture two point func- 
tion. Using the fact that the genus g contribution to the correlation function 

- y^ i ~ p ~ 1 k +l—p—2g 

depends on to according to < O ko ...O kp ^ 1 P > g ~ i=0 1 , we can 

write a formula where we exhibit the explicit genus dependence: 

-\{k - l)k{k + 1) < O k - 2 O kl ...O kp _ 1 P > 9 _! 
+4t (k + 1) < O k O kx ...O kv _ x P > a +At 2 {k + 2) < O k+2 O kl ...O kp _ 1 P > g 
+2 Ef=i(2A; + h + 2) < O k+ki O kl ...O kr -O kp ^P > 9 = 0, 

(7.12) 

where the first term is understood to be zero for g = 0. 
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These recursion relations are similar to the ones that exist for c < 1. 
The genus mixing first term looks like a contribution from a pinched handle, 
while the last term corresponds to contact terms between operators. What is 
lacking are the quadratic terms which arise when the surface is pinched into 
two separate surfaces. The reason these are absent is that we have considered 
1PI amplitudes with respect to the puncture which precisely excludes such 
terms. The corresponding recursion relations for connected amplitudes will 
be derived in the next section. 

It is straightforward to check that the results for the sphere, obtained in 
[p9[| using different methods: 




<O kQ O kl ...O kp _ x P> Q =-\hxii\[ (7.13) 



for even fcj's, satisfy ( |7.12|) , where E = |2f= ^i- Note that the first term 
of ( |7.12| ) does not contribute. The factors of t 2 assure that t 2 derivatives 
generate 2 insertions. Note that, in using the recursion relations, we have 
to fix the normalization < PP >o= — ^|ln/i| for t 2 = —2 (0/ = 1/4). 
The formula is also true for odd fcj's, unless 2E is odd, in which case the 
correlation functions are zero. This follows from (|7. 1 1|) for k — — 1, p = 1 



which states that < 0\P >= 0, as we discussed in the previous section. 

Let us explicitly evaluate a simple example of a correlation function at 
arbitrary genus. We consider the operator 2 . If we put k = in (|7.12| ) the 
first term vanishes, i.e. no genus mixing, and we are left with 

8t 2 < 2 O kl ...O kp _ 1 P > g 

+4t < O^.-.O^PP > g +(4(p - 1) + 4E) < O^.-.O^P > g = 0, 

(7.14) 

i.e. 

2t 2 < O^.-.O^P > g = -(p-l + E + to J-) < O kl ...O kv _ x P > 9 . (7.15) 

Since < O kl ...O kv _ L P > g ~ t^ +2 - p - 2g we find 

2t 2 < 2 O kl ...O kv _,P > g = -(2E + 1 - 2g) < O kl ...O kp _ x P > g . (7.16) 
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This expression is similar to the dilaton equation obtained in |26| for topo- 
logical gravity. In particular 

n 

< 0^02 > 9 = G^T f[(2 P - 1 - 2g) <> 9 . (7.17) 

^2 P =i 

Note that < O2 > g = ^(^g — 1) <> g , which suggests that 2 is almost, but 
not quite, the dilaton. 



7.4 Connected Recursion Relations 

Let us now derive recursion relations for connected amplitudes. This is not 
more difficult than in the 1PI case. In the case of the 1PI (with respect 
to the puncture) amplitudes the generating functional, T, depends on the 
nonpuncture sources t k (k > 1) and the nonamputated puncture one point 
function t = (3fi. t is kept constant while the other t k s vary. The generating 
functional for connected amplitudes, the vacuum energy E, depends instead 
on the t k (k > 1) and A, the puncture source, i.e. the cosmological constant. 
In this case A is kept constant when the t^s vary and to is a function of the 
tfc's. A is also the amputated one point function. T and E are related by the 
Legendre transform 

E(A) = (3 2 Ah-T(/j,). (7.18) 

We have 

dT dE d 2 E d/j d 2 T 

' ~Wiv ^~w^y dt k d(pA) ~~dKdt k d(P(j,y ( - } 

as a consequence of 

da du du dA 
n = — = — — -i — 7 20) 

dt k dt k 3Adt k K ' J 

The 1PI and connected two point functions are simply proportional to each 
other. The same recursion relation is valid. The difference comes when we 
take t k derivatives to obtain higher point functions. For connected diagrams 
to = /3fi is no longer to be kept constant. The resulting recursion relations 
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are no longer linear but rather quadratic 

-\{k + l)k(k - 1) < O k . 2 O kl ...O kp _ 1 P > c +2t 2 (2k + 4) < O k+2 O kl ...O kp _ 1 P > c 
+A{k + 1) £ iefc < O h ...O in P > c < O k O in+1 ...O iv _ x P > c 
+2 Ef=i(2A; + h + 2) < 0^ ki O kl ...O ki ...O kp _ 1 P > c = 0. 

(7.21) 

Note that £ =< P >c- We recognize an extra term suggestive of a surface 
pinching into two. If we extract the explicit genus dependence we find 

-\{k + l)k(k - 1) < O k ^O kl ...O kp _,P > Ci9 _! 
+A(k + l)E 9l+92=5 E i6fc < O h ...O in P > Ci9l < O fc O iti+1 ...O ip _ 1 P > Ci92 
+4t 2 (fc + 2) < O k + 2 O kl ...O kv _ x P > c , g 
+2E?=i(2A; + A* + 2) < O h + ki O ky ...d ki ...O kp _ 1 P > c , g = 0. 

(7.22) 

This structure is very similar to the one encountered in c < 1 theories [p6fl . 

All connected amplitudes may also be built by using the 1PI building 
blocks. Let us consider a simple example. Consider 

< O k P > c = -|£ < O fe P >= -- < O fc P > . (7.23) 
aA p 

The 1/p prefactor simply adds an external propagator as it should. Now take 
a ^ \A=const derivative. We get 

<o k o lP >c=(^f + ^|t)<O fc P> 

(7.24) 

-±(< O fc O,P > +§£ <O k PP>). 
Now use |f + f^f£ = to obtain 
< O fc O ; P > c = 4, < OiPP >< PO k >-^< O k P >< PPP >< POi > 

-i < O t O k P >+jz< OiP >< PPO k > . 

(7.25) 
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It is easily checked that this expression indeed obeys ( [7.2 1| ) using the recur- 
sion relations for the 1PI amplitudes. The 1/p factors provide external and 
internal propagators where needed. 



7.5 Virasoro Constraints 

There is an interesting way of rewriting the recursion relations ( |7.10| ), which 
explains where they really come from. They are simply equivalent to 

L k < P >= \/k > -1, (7.26) 

where 

1 8 °° d 

L k = —Ak + l)k(k - 1)^— + £ t p (2k + p + 2)-—. (7.27) 

4 Otk-2 p=Q (Jtp+k 

It is straightforward to verify that the L k s obey a Virasoro algebra 

[L n , L m ] = (n- m)L n+m . (7.28) 

The Virasoro generators can be understood as representing the algebra of 
coordinate transformations in the Liouville coordinate, i.e. ~ 
Let us see how this works. In the full quantum treatment, it is convenient to 
use the matrix eigenvalue A and its conjugate momentum p. The generator 
of the coordinate transformation is then pX k . One can convince oneself that 
one need not worry about different orderings in this case. Acting on the 
Hamiltonian we find two parts, 

8TI 

[p\ k ,U(X)]=X k — (7.29) 

and 

[pX k ,p 2 ] = k(X k -y + p 2 X k - 1 ) - ^(k - l)(jfe - 2)A fe " 3 . (7.30) 
When this is evaluated at the Fermi surface in accordance with the variational 

2 

prescription of Q4.23Q and ( |4.24j ), the ^- gives a U(X). Hence 



1 s-)T T 

< {—k(k - l)(jfe - 2)A fc - 3 + 2kX k - 1 U{X) + X k —)P >= (7.31) 
4 oX 
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which is precisely equal to ( |7.10| ). 

Unfortunately, the L_i constraint for the quadratic critical point is quite 
trivial 

L_! < P >= 2t 2 < dP >= 0. (7.32) 

It does not seem to yield a string equation determining < PP >, unlike the 
case for the c < 1 models, Iffifl . 

A natural generalization is obviously to consider general transformations 
p n X m . As far as X k correlation functions goes, this will not give anything new. 
The constraints one obtains are simply linear combinations of the old ones. 
Although, of course, they represent the full algebra not only the Virasoro 
part. 



7.6 Summary 

In this chapter we have seen how to obtain recursion relations for both 1PI 
and connected zero momentum correlation functions. With the help of these, 
any zero momentum correlation function for any genus may be calculated, 
once we know the puncture two point function. The recursion relations were 
furthermore found to correspond to coordinate transformations in the Liou- 
ville coordinate. This was the reason why our recursion relations could be 
expressed as constraints obeying a Virasoro algebra. The presence of the 
Virasoro algebra, and as we already have noted, a full Woo algebra, is of 
utmost importance. We will in chapter nine and ten take full advantage of 
this algebra. 

By adding time independent operators to the potential the critical behav- 
ior may change. The recursion relations and their Virasoro representation 
should be useful for studying such other critical points, different from the 
ordinary quadratic. They should also describe the flows between different 
critical points. Although it is not clear how to identify these new theories, 
it is easy to extract some information about them. The Virasoro constraints 
which we have derived are valid for any potential and provide relations be- 
tween different correlation functions. 

As we have seen, the recursion relations do not determine < PP >. 
This has to be specified independently, before the recursion relations can 
determine all other correlation functions. In general, when we consider some 
nonquadratic potential, more data need to be specified. For a theory with 
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a A 2 ™ as the leading term in the potential, we will need to specify < PP >, 
< PO\ >,...,< POm-i >• The correlation < P02 n -i > is always zero. The 
scaling of < PP > can be derived directly from the WKB approximation 
and is given by /i 1 / 2 ™ -1 / 2 . 

So much for zero momentum correlation functions. In the next chapter 
we will generalize to nonzero momentum. While doing so, the importance of 
the special states will become even clearer. 
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Chapter 8 

Nonzero Momentum 
Correlation Functions 



8.1 Introduction 

Let us now consider the matrix model and its special states, also for nonzero 
momentum. 

As we have seen the natural candidates for special state correlation func- 
tions are correlation functions of powers of the matrix model eigenvalues. 
Such objects were studied in |L5|] and the expected poles for discrete mo- 



menta were found. In this chapter we will extend some of the results of |15| 
and the previous chapter, for Ward identities. In particular, we will find that 
the recursion relation (|7.10| ) is in fact the Wheeler de Witt equation for the 
two dimensional world sheet quantum gravity. 

The recent developments revealing the symmetries indicate however 
that this is not the whole story. One should also consider correlation func- 
tions involving powers of the conjugate momentum. This we will do in the 
next chapter, where we will show the existence of this symmetry which 
will greatly simplify the subsequent calculations. 



8.2 The Wheeler de Witt Equation 

Ward identities for correlation functions are in general obtained by changes 
of variables in the path integral. Examples of such Ward identities were 



62 



obtained in the last chapter from simple coordinate changes in the matrix 
eigenvalue. They can be thought of as generated by commutators, or classi- 
cally, i.e. on the sphere, Poisson brackets with pX m . They obviously obey a 
Virasoro algebra, which is part of a algebra generated by all monomials 
p n X m . We may also introduce time dependence and consider generators with 
certain momenta q, i.e. p n \ m e % it _ L e ^ us, as an example, make a pX k e iqt 
variation of the one puncture function. 

As we saw in chapter three, the two puncture function is schematically 
given by 

<PP>=-lm [°° dT [ dXG(X,X;T), (8.1) 
n Jo J 

where the calculation is done at the Fermi surface. We will shift its energy 
to zero, hence putting the Fermi energy as a constant term in the potential. 
G is the path integral given by 

G(A 1 ,A 2 )= f X2 [dpdX]e^I dt ^~ l 2P 2 +uW) y ( 8 . 2 ) 

J Ai 

in Euclidean time, where U(X) = J2 P t p X p is the potential. The variation of 
the partition function would have involved a sum over all states in the Fermi 
sea up to the Fermi level. By inserting a puncture, i.e. taking a derivative 
with respect to the Fermi energy as in the last chapter, we restrict ourselves 
to the Fermi surface. Next we perform the variation of < P >. The measure, 
as given by (|8.2| ), is invariant under the change of variables. (Clearly, we are 
not supposed to differentiate with respect to t when changing variables in 
the measure.) The change in the action gives rise to the following identity 
among two point functions 

< ( / dt(iqX k p + kX k ' V + ^rX k )e iqt )T >= 0. (8.3) 
J oX 

The puncture is now a tachyon, T, carrying away the momentum. The piece 
with a single p is evaluated by integrating over p in the path integral, obtain- 
ing a A which then is partially integrated. We then switch to a Hamiltonian 
formulation, remembering that we should use Weyl ordering. We finally 
obtain 

-h(k-l)(k-2)<O k - 3 P> g -i +Y,( 2k + P) t p < P+ k-iP > g = 0, (8.4) 

p 
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to all genus at zero momentum, as in the last chapter, and 



a 2 

£(2A; + p)t p < O p+k ^T > g>g < O k+1 T > 9i5 = (8.5) 

on the sphere for nonzero momentum. The combined case is somewhat more 
complicated due to the need of taking an imaginary part and an effective 
imaginary shift iq in to, relevant for higher genus. We will not consider this 
any further. In these equations q indicates momentum and g genus. We have 
introduced the notation O k for X k . The different genus for the first term in 
( |8.4p is due to a p, A commutator which arises when we want to evaluate the 
p 2 against the wave function of the Fermi surface. This gives an h, which is 
the same as the genus coupling constant. Following we may define the 
loop operator given by 

w(l) = e lx . (8.6) 

It corresponds to cutting out a hole in the surface with a boundary of length 
I. The reason is the following. If we insert a power n of the original matrix 
variable <fi on the surface this creates a little hole, the length of the boundary 
being proportional to n (the number of legs) and the lattice spacing a. For 
fixed n the length clearly shrinks to zero in the double scaling limit. To get a 
finite length we must also take n to infinity. Introducing A as the eigenvalues 
m of expanded around the top of the potential, we find in the double scaling 
limit 

m n ~ (1 + a\) l/a -> e lx . (8.7) 

We may then Fourier transform to obtain a differential equation in the loop 
length. We get 

E M V + 2 l d^ ] + to/2 - /4] < w(l)T >= (8 - 8) 

at zero momentum and 

E ^ l2 w + + toP + q2] < w{l)T >= (8 - 9) 

on the sphere. The / 4 term is of order one higher in the string coupling 
and do not contribute on the sphere, to is the Fermi energy fi with the 
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appropriate number of /3's absorbed. In the case of the usual harmonic 
oscillator potential, where to = the resulting equations are in fact the 



Wheeler de Witt equations obtained in [58]. There they were inferred from 
the explicitly computed two point function. 

On the sphere this is one of the most striking verifications of the equiv- 
alence of the Liouville and matrix models. At zero momentum, the Wheeler 
de Witt equation is just the Fourier transform of the Gelfand-Diiki equation 
for the resolvent for the Schrodinger operator. This was the way in which the 
zero momentum version of (|8.4| ) was derived in the last chapter. If we want to 



be careful we need to rescale A by y/—t2 to get a dimensionless I, see section 
9.4. This is needed for the exact correspondence between the above result 
and the mini super space canonical quantization of Liouville theory. Recall 
that, |40| , t 2 = — 2^7. From ( |8.9| ) one might try to draw some conclusions 



about the Liouville theory correspondence to the more general potentials 
above. Clearly the last term, which corresponds to the matter piece, does 
not change while we change potential. Instead, it is the piece which would 
be expected to arise from a canonically quantized kinetical term for the Li- 
ouville mode, which gets modified. Hence one is led to the conclusion that 
these more general models (however with p independent potentials) may cor- 
respond, not to modifications of the matter theory, but rather to different 
theories for the Liouville part. This is also consistent with the point of view 
for which this paper will argue, that the special states must be represented 
using both A's and p's, not just the A's. 



8.3 Two Point Function on the Sphere 



So far we have seen how to derive some simple recursion relations for various 
correlation functions of the matrix eigenvalue. We will now turn to some 
other methods for calculating correlation functions. We will obtain results 
for all genus, but we will begin by deriving an expression for the two point 
function of special operators with nonzero momentum on the sphere. 

We will use the same method as in 
given by 

1 



40 . The two point function is then 



< OrOe > 



(2^) 2 J? 



n 2 uj 3 



j p 2 + n 2 uj 2 



cos nujTi cos nuT 2 , 



(8.10) 
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where T = c l\[ot'\ log /x| is the period of the classical motion and uj = 2n/T. 
t is the classical time of flight variable, which for the present purposes may 
be though of as the Liouville mode. We have chosen r to be zero at the 
classical turning point. This is the natural choice since the critical behavior 
arises from this region. The classical equation of motion at the Fermi surface 
is given by 

^ = fifo - U(X) hrU(X) = ^(A 2 - A 4 ). (8.11) 

If we expand about the top of the potential, A = x + keeping only the 
quadratic term, we get 

= V2rtfi 1/2 cosh -J=. (8.12) 



Hence 

O r = x r = (2a') r/ V /2 cosh r -4=- (8.13) 

Va' 

If we now tried to use (|8.13|) in (|3.10|) and perform the Fourier transforms 




we would get in trouble. As /i — > and T — > oo, (|8.13 ) diverges at the limits 



of the integration. The reason is that our approximation is good only up 
to \t\ = T/4 and breaks down thereafter. One way to take care of this is 
to simply put in a cutoff. The critical behavior should not depend on the 
nature of the cutoff. Let us however be slightly more careful and study the 
exact solution of (|8.11|) in terms of elliptic functions. With our convention 
for r, the solution is 

A(r) = A_crf(^^r) A| = \ ± 2v / ^ 7 /i 1/2 . (8.14) 

\/2a' 2 



The function cd is the same as sn translated a quarter of a period, and can 

cn(u) 
dn(u) ' 



be written as cd{u) = T^l , where for small \x the elliptic functions cn and 



dn are given by |l 



and 



1 mi 2 u sinhw 

cn(u) -2 -(tanh u -3—), (8.15) 

cosh u 4 cosh u 



, . . 1 mi . , 9 u sinh u . . 

dn(u) -^ + — (tanh 2 M + -5-), (8.16) 

cosh u 4 cosh u 
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where m\ = \ — (t— ) 2 . 



• i ji ■ i / \ 1— —r- (sinh 2 u— u tanh u) 

Consider the expansion cd{u) ~ 1+ A (sinh . M+MtanM , 



valid to only first order in m\. For ^(sinh w+wtanh-u) < 1, we may expand 
the denominator to get 



cd(u) ~ 1 H — — !- cosh 2m. 

4 4 



(8.17) 



From this we recover (|8.12| ), but now we see that it has limited validity. As 
fi and mi go to zero, the expansion can not be trusted above |r| = T/4. 
Performing the Fourier transforms up to |r| = T/4 we obtain 



= (r/2 - fc) sinh(r/2 - k) 



2VV 



fc=0 



7 (r/2- k) 2 + n 2 u 2 



-(-l) n/2 , (8.18) 



for n even and 



,a 
"~2 



rV /2 E 

fc=0 



r 
k 



2nu cosh (r/2 — k) 



2vV 



±(r/2-k) 2 + n 2 u< 



l(_l)(n-l)/2 ) ( 8 _ 19 ) 



for n odd. To compare with previous results, we assume the same contri- 
bution from both tops of the potential (around r = and r = T/2). This 
doubles the Fourier transforms. We then need to rewrite the sum over n as 
a contour integral, i.e. 



1 a \(r+s)/2 n (r+s)/2 / dz Az 2 ^ ^ 



4tt 2V 2 



x cot 



7T2 



4^ 2 + z 2 fc=0i=0 v 
16(r/2 - k)(s/2 - I) 



r 
k 



2oo (4(r/2 - k) 2 + a'z 2 )(A(s/2 - I) 2 + a'z 2 ) 



x sinh 



r T 

- - k) 

'2 >2^fti 



— tan 



7T2; 



sinh 

4a'^ 2 



2 >2sfd 



2u (4(r/2 - k) 2 + c^ 2 )(4(s/2 - I) 2 + a'^ 2 ) 



x cosh 



r --k) — 
2 J 2Vrt 



cosh 



2 ^v 7 ^ 7 



(8.20) 
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where the contour wraps around the real axis. This is obtained by using 



£ H{n) = - f dzcot(^)H(z) £ H(n) = -- f dzt a n(-)H(z), 

n=even 4Z J Z n=odd ^ l J Z 

(8.21) 

and a rescaling z — > z/u. 

We can rotate the contour in (|8.20|) to wrap around the imaginary axis 
and evaluate the expression by picking out the poles. The z = ±ip poles 
correspond to tachyon contributions while the other poles are associated 
with the special operators we are concerned with. There are both simple and 
double poles of this kind. The simple poles do not give rise to any singularity 
in fi. From the double poles we get | ln/x| terms from derivatives of the cot 
and tan factors. We finally obtain 

,a\ (•■+*) y/c/ (r+s) \ 4(2: _M2 

(_)._„, M T,{ k )[^_ k ) a y_ 2 4(i V (S-^) 

where we have assumed r < s. Note that the k = r/2 term only contributes 
at p = 0, where the above expression indeed agrees with previous results for 
the zero momentum special operators. Indeed, when p = we find 

-(^^^IM( (r r + + ^). (8-23) 

in agreement with (|7.13[) for t 2 = — ^7. We will return later to a discussion 
of ( |8.22| ) and in particular the significance of the poles. 



One might note that the rotation of contours is precisely related to the 
analytical continuation to a right side up oscillator which will be the basis 
for the calculation in the next section. 



8.4 Two Point Function for all Genus 

We will now derive a formula for the special operator two point function 
which is valid for all genus. The starting point will be the puncture two 
point function 

1 00 1 

< pp >= — im E ~p — r> ( 8 - 24 ) 

7T n=0 tin ~ t 
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where E n are the energy eigenvalues of the upside down potential and hence 
imaginary. We will be especially interested in the quadratic case where E n = 
^=^(2n + l)i and to — ftf^- However, the following derivation is valid for 
the general case. We will conclude this section with a few comments on the 
changes that occur when we have an arbitrary potential. 

Correlation functions of arbitrary, in general time dependent, operators 
are obtained by adding these to the Hamiltonian and evaluating the energy 
eigenvalues E n to the appropriate order in perturbation theory. This is a 
well defined prescription even for time dependent perturbations. To do this, 
we evaluate the Euclidean matrix element 

<n\T(e~f-^ m )\n> , (8.25) 

where H(t) = H + V(t) with H Q is the unperturbed inverted harmonic 
oscillator Hamiltonian and V(t) is some time dependent perturbation which 
we will take as €ie ipit Ok + e 2 e %V2t Oi. If we expand to second order in ei and 
e 2 we get 

t\t 2 / dt 2 / dti < n\Ok\m >< m\Oi\n > 

"2 

To do the integrals we extract a common factor e~ EnT and recall that the 
.En's are imaginary. Dropping oscillating terms, we find in the limit T — > oo 
(for momentum conservation pi + p 2 = 0) 

eie 2 e-^ T - 'L—. (8.27) 

E m - E n + ipi 



Since we have two different time orderings we obtain for (|8.25|) 

2(E/ T ra E" n ) 



e -E n T 



1 + exe 2 T < n\Ok\m >< m\Oi\n > 



m=0 



p 2 + [E m - E n ) 



(8.28) 



Interpreting this as e~ EnT (l + eie 2 SE n T) ~ e -{E n -eie 2 SE n )T ^ we g n( j ^ e per- 
turbation in E n . 

The above result then allows us to evaluate the second order variation of 
the two puncture correlation function to be 
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i Q2 oo i 

< PPO k O { >= ———Tm^- -| ei=£2=0 

7T de 1 de 2 ^ E n - t 

= ^Et/t^ < n\O h \m >< mmn > - E j (8.29) 

n n,m {En ~ *oJ V + \E m ~ -&nj 

or after removing one puncture, 
^ P n n ^ X t ^ < n\O k \m >< m\Oi\n > {E n -E m f 

< POk ° l >_ ~ Im i (K-t )(^-to) ? + {E m -E n r (8 - 30) 

(Recall that the En's are imaginary.) Let us now limit ourselves to the 
quadratic critical point, i.e. the inverted harmonic oscillator, and compare 
with our previous results for the sphere. 

The matrix elements are polynomials in n. On the sphere, only the highest 
power of n survives since this term will give the highest power of fi/3 after 
the sum is performed. Using < n|A|ra >= (^^f-) 1 ^ 2 (S nlm+ l^/n + 5 n+ ^ m y/m) 
the leading terms of all the matrix elements form a Pascal triangle, hence 

< n\O r \n + k >~ ( (r _\ )/2 ) n^C-^fr\ (8.31) 

for k = — r, — r + 2, r, otherwise zero. In the spherical limit we can replace 
the sum over n by an integral. We get 



[m) (r+s)/2 



Im — i — - r — ~ Re 

5^7(271+1)1-^ 



/-oo (j„\(r+s)/2 
/ g(r+.)/2 (v ^7)(r+,)/2+ 1 f 



x + m 

(8.32) 

The integral is defined by taking a suitable number of derivatives and keeping 
only the part singular in \i. With even r + s we get 

( M /3) (r+s)/2 (v / ^ 7 ) (T " +s)/2+1 | ln^|. (8.33) 

Odd r + s give zero, confirming that odd-even correlators are zero. Putting 
all of this together we recover Q8.22 ). 

As we have indicated ( |8.30| ) is valid for a general potential. The poles 
have a simple explanation. They appear at the values of momenta that can 
excite transitions between the energy levels in the upside down potential. 
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In the case of a quadratic maximum these are linearly spaced and lead to 
integer values (in units of of the momentum. If we, however, consider a 
critical theory generated by a potential whose leading term at the maximum 
is higher than quadratic, the situation will be more complicated. The energy 
levels in the upside down potential will now no longer be linearly spaced. 
The resonant values of the momenta will be differences of these levels, and in 
general we would expect a strange transcendental set of discrete momenta. 
What conformal field theory could produce such special physical states? 



8.5 Correlation Functions from Path Integrals 

In ||58|| , a generating functional for all correlation functions was calculated 
using path integrals. In this section we will indicate how our results for the 
special operators agree with this approach. We will consider both the sphere 
and later all genus. 

Let us first state the results of [Q. We write them as 

ar+s r oo ij3n£+i/2cath(Z/2)(zf+z%) 

<O r O s P>= - - Im/ de- 
o r z 1 o s z 2 Jo 



sinh(£/2) 



x 



+ E 

fc=0 



2iTe 



Ai k k 



T|p|/2 sinh(b|e/2) j 
sin 7r |p| 



Z\Z 2 



J, 



Z\Z 2 



k 2 -p 2 KV sinh(£/2) 



sinh(£/2)' 
) sinh(fc£/2) 



(8.34) 



-I 21=^2=0 

The derivation of this expression is very straightforward. One simply needs 
the propagator of the harmonic oscillator. One then need to Fourier trans- 
form to momentum space and also make an integration against the loop 
operators, 



3S 



The term which will be of interest to us is the last one. It corresponds to 
the contribution from the special states. Let us first consider the spherical 
case. We expand (|8.34j) for small £ and pick out the term of order q = (r+s)/2 
in z\ and z 2 . A small £ expansion will clearly give the string perturbation 
series. Doing so we obtain 



Im 



^z 2 )H« I di- 
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1 



2k—nn—2k 
Z-| z 2 



Al 2 



p 2 -I 2 



(8.35) 



The integral is defined by taking a suitable number of derivatives and ne- 
glecting terms analytic in /i. For fixed powers of both Z\ and z 2 we get 



z\z 2 r + s 
= r!s! 2 



(/3 At )( r+s )/ 2 - 1 |ln/x| 



fc=0 



r 
k 



4(r/2 - fc) 2 



(s + r)/2 - fc y p 2 - 4(r/2 - A;) 2 ' 



(8.36) 



Hence we find agreement with our results for the sphere apart from an unim- 
portant normalization factor. Let us also check the formula for general genus 
in a few cases. For simplicity we will limit ourselves to r = s. In that case 
we get 



Im 



(^)f e^sinh(M). Ak 



2 p 2 - k 2 



r — k 
2 



1 



r — I 



X S^-0! V (r~k)/2-l 



r — I 



r+l-2l 



{r + k)/2-l J v sinh£/2' 



(8.37) 

We may identify the harmonic oscillator matrix elements by expanding the 
inverse powers of sinh£/2 in powers of e~ n ^ using 

l sinhe/2 ; to 1 ' 

For k = r, which is the term with the highest momentum pole, we find 



Im 



,IZiZ 2 , 







°°^ s inh(^ — 



1 



\r+l 



2 p 2 — r 2 sinh^/2' 



5.39) 



To relate this to our previous formulae we use the recursion relation for Her- 
mite polynomials H n+ i(X) = 2xH n {\) — 2niJ n _i(A), and also the expression 
for the normalized wave functions i]) n (X) = 7r i/4( ra i)i/2 2 ^/2 H n (X)e^2 X . From 
which it follows that 



< n\X r \n + r >= ( ( n + r )"-( n + 1 ) )i/2 ( _ tt > ) r/4 



5.40) 
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We then obtain (with a' = 1) 



2 r+3 Im 



riri Jo 



1 - e~ r? ) < n\X r \n + r > 2 e" ( ' 



2 -( n +l/2)£ r _ 



n=0 



i) 



Similarly, for k = r — 2 we find 



2 r+3 Im 



jPrt sinh 



(r-2)A r-2 



p 2 — (r — 2) 2 



" 2 , 1 



x 



r! sinh | 



j (r-l)! l sinh| ; , 



(8.42) 



Now using 



< n | A > + r - 2 >= ( (n + r 2 )-( w + ^ jV^ + L_±)(- a >y/* ( 8 . 43 ) 



which again is a result of the Hermite recursion relation, we can rewrite ( |3.42j ) 

as 

T ! 



2 r+3 Im 



1^1! r e <^(l_ e -(r-2)^ 

~!r! Jo 



X 



]T <n|A> + r-2> 2 e~ (n+1/2)5 : 



r-2 



8.44) 



P 2 " (r - 2) 

To compare with the previous expression, we rewrite our formula ( |S.30| ) as 



— Im 



;^(e~ - 1) 



x 2^ < n|Oi|n + >< n + fc|0 2 |n > e s 



n=0 a 'P 2 - k2 

which indeed, with a' = 1, agrees with the examples above. 



.45) 



8.6 Identification of States 



In this section we will compare the matrix model results with the predictions 
of the continuum theory. We will try to identify the matrix model operators 
which can be used to excite the continuum theory special states. 
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Consider a correlation function < O n O n > in the matrix model. Accord- 
ing to ( |S.22| ) the scaling dimension of the O n operator is d = |. Further- 
more, we have noted the presence of poles at momenta p = ±7^7=7? where 
k = —n, —n + 2, ...,n. These combinations of scaling dimensions and mo- 
menta are the same as those which we find in the continuum theory. As was 
explained in the introduction, there exist in the continuum theory operators 
with dimension d = and momenta p = corresponding to level rs 



2 

descendants of tachyons at that momentum. We are therefore tempted to 
identify the p = ±^7=7 pole of < O n O n > as coming from the level \{n 2 — k 2 ) 
descendant of the p = ±77=7 tachyon. It is reasonable to assume that the 
matrix model operators O n with these momenta can be used to excite the 
continuum special states according to the above identification. The zero mo- 
mentum states, however, differ from the above scheme. There are no poles at 
zero momentum. On the other hand, it is only for zero momentum which the 
term in ( |3.30| ) with diagonal harmonic oscillator matrix elements contribute. 

Let us illustrate the above procedure by a couple of examples. The 2 
operator has scaling dimension d = 1 and the < O2O2 > correlation function 
a pole at p = ±7^=7- in the continuum theory this combination of dimension 

and momentum corresponds to the special tachyon operator e ±l ^ fX . A more 
interesting example is O4. It has scaling dimensions d = 2 and its two point 
correlation function poles at p = ±7^7=7 and p = ± 7^7=7 • The pole of highest 
momentum is again one of the tachyons with descending null states, this 
time e ^ , but the other one would seem to be related to the level 3 

descendant of momentum p = ±7^=7- We can also consider odd operators. O3 
has scaling d = 3/2 and, in addition to a tachyon pole, a pole at p = ±-^7=7, 
which can be identified as the level 2 special state given by (8). 

We may also consider correlation functions between different matrix model 
operators, i.e. < O n O m > with n^m. The number of poles is now deter- 
mined by the lowest dimension operator. Since we have nonzero correlation 
functions between different operators there is substantial mixing. This means 
that the identification of matrix model operators and continuum theory spe- 
cial states is not one to one. This problem will be resolved in the next section 
where we introduce the matrix model generators of the symmetry. 
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8.7 Summary 



In this and the previous chapter we have studied the special operator correla- 
tion functions of c = 1 quantum gravity. For nonzero momentum, the matrix 
model correlation functions of time dependent perturbations of the potential 
about its maximum are found to have poles at the special values of momenta 
predicted by Liouville theory. The operators also have the correct scaling 
dimensions. This allow us to identify the matrix model operators which may 
be used to excite any of the special states in Liouville theory. 

We have also derived expressions for the two point correlation functions 
at nonzero momentum. The correlation functions which we have derived are 
found to have poles which, for a general time independent potential, are given 
by energy level differences in an anharmonic oscillator. In the harmonic case 
the levels are equally spaced and precisely coincide with the special momenta 
of Liouville theory. In the general case the meaning of our results are not 
clear. 

We could also imagine adding time dependent pieces to the potential. 
Clearly, this will drastically change the physics if the time dependence is such 
that we sit on top of one of the momentum poles. We have, in chapter 4, 
speculated that this leads to models in which the target space is punctuated- 
i.e. the continuous line is transformed into a set of discrete points with 
spacing related to the special values of momenta. It would be very interesting 
to explore the structure of these theories. 
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Chapter 9 

The Matrix Model W { 
Symmetry 



9.1 Introduction 

In previous chapters we have seen how to compute certain correlation func- 
tions of the matrix model eigenvalue. We have done both direct calculations 
and derived recursion relations. There is nothing wrong with deriving re- 
cursion relations in the way we have done. In particular, one can certainly 
extend the results to correlation functions involving powers of the conjugate 
momentum. However, it is more convenient to make use of the large set of 
symmetries in the theory, the W^. In this chapter we will see how this is 
realized in the matrix model. 

After we have done that we will use the algebra to calculate correlation 
functions. We will also compare with the results obtained in previous chap- 
ters. 

The calculations will focus on the sphere. Higher genus are reserved for 
the next chapter. 



9.2 The Woo Algebra 



As shown by Witten in [73[] it is convenient to change basis to (A— p) n (X+p) m 
(at a' = 1) rather than dealing with the A and p separately. For certain time 
dependent prefactors, we then obtain transformations which leave the action 
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invariant. These transformations are generated by, in Minkowski time, 



W r ' s = (\ + p) r \X-p) s e qt (9.1) 

with q = r — s. These obey the algebra 

{W ri > Sl , W r2 > S2 } = 2(r 2Sl - r lS2 )W ri+r2 - l ^ +S2 -\ (9.2) 

classically generated by the Poisson brackets {X,p} = 1. For a general 
momentum q in ( |9.1| ) we find, when acting on the Minkowski action S = 
/(p A-i(p 2 -A 2 )), 

{W r ' s ,S} = (r-s-q)W r > s . (9.3) 

Hence a symmetry for appropriate discrete values of imaginary momentum. 
This is equivalent to saying that W = W(p, A, t) is a solution of 

dW dW 

_ = W + {H , W] = . (9 . 4) 

Expressed in terms of the initial conditions, po and Ao, we have W = W(po, Ao) 
i.e. any time independent function of the initial conditions. The generators 
Q9.1D are then simply obtained through evolution in time. Hence the trans- 
formations can be understood as time independent canonical transformations 
of the initial conditions. 

A more convenient way of labeling the generators is through their SU(2) 
quantum numbers J = (r + s)/2 and m = (r — s)/2. With the definition 

Wj, m = (A + p) J+m (\ - p) J ~ m e 2m \ (9.5) 

one gets 

{W Jumi , Wj 2 , m2 } = 4(m 2 Ji - m 1 J 2 )W Jl+ j 2 

-1,1711+1112- 

(9.6) 

As in the previous chapter, we will continue to Euclidean time, i.e. t — > it, 
and consider the system as a right side up oscillator continued to imaginary 
frequency. The latter effectively means p — > ip and an extra % in the structure 
constant of ( |9.6| ) and the eigenvalue of (|9.3j). As already emphasized several 
times, this analytical continuation has several advantages. Contrary to the 
case above the elements now act like step operators capable of exciting 
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Figure 9.1: J = 3/2 and J = 2 special tachyons. 



the system. The operators (|9.5|) however, make a step in imaginary energy 
and in fact take you out of the Hilbert space. This has also been discussed 
in 



19 



The continuation also leads to the possibility of drawing pictures of the 
special states. This was mentioned at the end of chapter 4. Recall figure 4.3. 
With the identification of the special states above as generators, let us 
consider a couple of examples. The actual special state perturbations would 
look like 

(p 2 + X 2 ) J - m ((X + ip) J+m + (A - ip) J+m ). (9.7) 

Both signs of the momentum need to be included. In figure 9.1 we find 
pictures of the J = 3/2 and J = 2 special tachyons. 



9.3 Correlation Functions from 

In this section we will calculate correlation functions involving the operators 
Woo as defined in the previous section. The symmetries will help us 
organize the Ward identities. As an example, let us start with the two point 
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function. Using 

1 00 1 

<PP>=—l m J2- (9.8) 

7T n=0 tin ~ t 

and simple perturbation theory we get 

„, „ r „ 1 T ^ ^ < n I Wi I k >< k I W 2 I n > 1 
< WiW 2 P >= Im 2^ 2^ 



71 timeord n,k E n ~ *0 iPl + ^fe ~ E n 

(9.9) 

Since the W's are of the form ( |9.5|) , continued to Euclidean time, they are 
simply raising or lowering operators in the continued harmonic oscillator. 
This means that only a few of the matrix elements are nonzero. Since W 
raises by 2m = r — s we get 

< w lW2 p >= _i faE <HM|n> 1 

7T V E n -t z(p 1 -2m 1 ) 

We have reduced the two point function to a one point function using the 
commutation relations. If we restrict ourselves to the sphere, use the algebra 
given by ( |9.2|) and directly calculate the one point function we get 



< Wl W 2 >= _ Mrn,Jr- rn,J 2 ) l^ 2 „ + „_, 
Pi - 2mi 7r Ji + J 2 



or equivalent ly 

— ?m, 1 

< >= (2/i) Jl+j2 I log// I . (9.12) 

2m! — Pl 7T 

The simplest way to obtain the one point function on the sphere is to use 
the classical Fermi liquid picture introduced by Polchinski H67|| . As we saw 
in chapter three, we simply need to do the phase space integral 

<W>= J dpd\W(p,X). (9.13) 

This was also discussed in [19 . To make everything well defined, we need 
to introduce an extra puncture, i.e. take a derivative with respect to the 
cosmological constant. Doing that, the integral over the whole Fermi sea 
becomes just an integral over the Fermi surface 

< WP >= <fw(p,X). (9.14) 
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The integral is to be performed along a hyperbola A ~ p = //. Although the 
answer is really infinite, we know that the piece nonanalytic in // is - | log \i \ 
for < PP >, i.e. when we only integrate 1. The other zero momentum 

\ 2 2 

operators simply involve integrations of ( — rp-) ra = /i n , again constants along 
the Fermi surface hyperbola, so we get 

< W n ' n P >= --{2n) n | log// | (9.15) 

7T 

and from this ( |9.11| ) follows. Our conventions are such that a' = 1. An- 
other approach, which is convenient when calculating correlation functions 
of nonzero momentum, is to continue to the right side up oscillator where the 
Fermi surface is a circle. In that case, however, we need to remember to put 
the Liouville volume | log/z | in by hand. Parenthetically we may note how 
a general correlation function may be obtained in this way. For instance, 
the two point function is obtained by perturbing the Hamiltonian, and hence 
the Fermi surface by one of the operators. If we integrate the other operator 
against the change in Fermi surface we get the correlation. 

Let us now consider the more complicated case of a three point function. 
In this case there are six different time orderings. Again perturbation theory 
gives us 

1 T x , ^ < n\Wi\ m >< m\W 2 \k >< k I W 3 I n > 
— 1m 

timeord n,m,k n ^0 

x . (9.16) 

ip x + E m - E n ipz + E n - E k 

Let us make the sum over time orderings more explicit. We find three terms 
of the form 

1 < n | (WiWWs + W 3 WWi) I n > 1 

— 1m } { — -+perm\. 

vr V E n -h i{pi - 2mi)i{p 3 - 2m 3 ) 

(9.17) 

This may, after some straightforward manipulations, be rewritten as 

< Wl w 2 w 3 >= <m&WM> + <W*WlW,]1> 

(pi - 2m l )(p 3 + 2mi + 2m 2 ) (p 2 - 2m 2 )(p 3 + 2m a + 2m 2 ) 

(9.18) 
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For the sphere, we now use the algebra ( |9.6| ) and the explicitly calculated 
one point function to get 

< W X W 2 W 3 >= 
(pi - 2mi)m 2 m 3 J x + (p 2 - 2m 2 )m 1 m 3 J 2 + (p3 - 2m 3 )m 1 m 2 J 3 
(pi - 2m 1 )(p 2 - 2m 2 )(p 3 - 2m 3 ) 

x -(2/i)^- 7 "- 1 I log// I . (9.19) 

7T 

The general higher point function can be obtained recursively from the 



three point by use of fl9.2|) and (|9T3| ). To get the N point function with an 
additional operator Wn, we vary the N — 1 point function with Wn knowing 
that the total variation is zero. The variation consists of two terms. One 
from varying the action as given by ( |9.3| ) and a sum of terms from varying the 
other operators as given by (|9.6| ). Each of the terms in this sum is obtained 
by shifting pi — > pi + p^, rrii — > rrii + and Jj — > Jj + J^+i — 1. We also 
need to multiply with the Clebsch-Gordan coefficient 2(Jjmjv — J^rrii). It is 
an easy exercise to check that (25) is obtained by applying this procedure to 
(|9~TlD or Q9T2D . 



There is one subtlety in the variational procedure which should be noted. 
Previously we derived recursion relations for zero momentum correlation 
functions by varying < P >. Recursion relations for < PP > are obtained 
from these just by taking a to derivative. Since the < P > recursion relations 
involve an explicit t we find an extra term, 2k < Ok-\P >, in the < PP > 
recursion relations compared to naive expectations. The reason is that we 
insisted on rewriting any conjugated momentum p in terms of A's. The rela- 
tion involves the potential U and hence depends on the coupling constants, 

2 

in particular t . The statement is that although < U(X)P >=< yP >, we 

2 

have in fact < U(X)PP >^< \PP >■ In the recursion relations we are 
discussing, we never perform any evaluations like this and hence need not to 
worry about these issues. For instance, if we want to include extra punctures 
in the correlation functions we are free to do so without extra terms in the 
Woo relations. 

Using this method one can write down several different recursion relations. 
One simple example is: 

< Tj yJ Wj umi Y[T JhJi >= — ^ \ < Wj+ji-ij+m! n T JiyJi > . 

i=2 P ZJ i=2 

(9.20) 
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We will come back to this relation later, when we compare with the Liouville 
model results. 

Rather than considering these general expressions, let us look at a couple 
of important examples where the form of the general iV point function is 
particularly simple. 

The first example is the N point function of special tachyons. It is given 

by 

N _ y-rN 9 t 9 ,jN-3 w 



71=1 



nt 2 (2^n-p„)vr^- 3 



The quantum numbers have been chosen as m n = J n for n > 1 and mi = 
— J\. This is just the pole part of the general tachyon correlation function 



as computed both in the matrix model |36|, [58) and in the Liouville theory 
[j36l , |24j| , up to a factorized normalization factor. The proof is by varying the 
three point. We cannot just vary the three point tachyon correlation function, 
since some of the J's are really m's in disguise and J and m vary differently. 
Instead we start with the general three point and make an arbitrary number 
of tachyon variations. A simplification is that we at each step only have 
to vary the single negative chirality tachyon. It is only from there we will 
get a nonzero Clebsch-Gordan coefficient. Following the prescription above, 
performing N — 3 variations we get 



N 



< n T n >■- 



N N 

(pi - 2m 1 + J2(Pi ~ 2m ! )) m 2%( h + J i ~ N + 3 ) 

i=4 i=4 



n=l 

N N 

+ (p 2 - 2m 2 ){m 1 + ^2mi)m 3 J 2 + (p 3 - 2m 3 )(mi + ^m i )m 2 J 3 

i=4 i=4 

x Uf=4 4( JiK + ZjLj+i mj) - ( Jj + Ef =w J,-N + i) mi ) 
(px - 2m! + E,= 4 fe - 2m,)) nt 2 (Pn - 2m n ) 

x -(2/i)Sn=i^- Ar + 2 | log// | . (9.22) 

7T 

The product in the denominator is the product of all the Clebsch-Gordan 
coefficients of the variations. Note that each get shifted by the successive 
variations. By the use of momentum conservation and evaluating the m's as 
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Figure 9.2: Graph for four point tachyon correlation function. 



J's, the formula (|9.21 ) is proved. This derivation shows how the combinato- 



rial factor from the /i derivatives is a consequence of the symmetry. 

Since this derivation is quite illuminating, it might be reasonable to con- 
sider a particular example which is more transparent than the general deriva- 
tion above. Let us consider a four point tachyon correlation function. A 
sketch of the relevant setup with three vertices and a one point is given in 
figure 9.2. The vertices give rise to the factor 

(J 2 m! - Jim 2 )(J 3 (m 1 + m 2 ) - (Ji + J 2 - l)m 3 ) 

x (JiiiTii +m 2 + m 3 ) - (Ji + J 2 + J 3 — 2)777,4) 
= 2Ji J 2 J 3 J 4 (Ji + J2 + J3 + J4- l)(Ji + J2 + J3 + Ja - 2). (9.23) 

The last factor is then cancelled by the one point function. The remainder, 
apart from a factorized part, is the well known Ji — 1 found by many other 
methods. 

If we want to consider the zero momentum operators, we have to be 
careful. The Clebsch-Gordan coefficients are zero in this case but these zeroes 
precisely cancel the momentum poles and leave a finite result. We get 

< LI W n >= — ^Zi(2/i) E ^ Jn I log// I • (9-24) 

n=l 71 a ^ 
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We will use induction for the proof. We find 



N+l j 

< IT W n >= -J2^(JN+im k - J k m N+1 ) 

n=1 k PN+l — ^N+l 

1 d N ~ 2 



X nd^^^^ 1 1 bg/i 1 • (9 ' 25) 

If we then put Pn+i = and use that the sum of all m's must be zero the 
result follows. This can also be checked by an explicit phase space calculation. 



9.4 Some Consistency Checks 

Let us now study the results of the previous section more closely. In particular 
we will compare the two point correlation functions with the results derived 
in the last chapter. These were correlation functions of pure powers of the 
matrix eigenvalue, without any momentum powers. In terms of the Ws they 
are given by 

On = t ( I ) W^ fe -"i (9.26) 
From this follows that the two point function is given by 

< O n O m > q = JL ± f ( I ) ( 7 ) < Wr k ' k W- 1 ' 1 > . (9.27) 

Using ( |9.12| ) and some simple algebra we find 

-1 1 i (w+m) . . A (n \( m \ 4(f-fc) 2 



In precise agreement with [I5| and the result ( 8.22| ) of the previous chap 



ter, recalling our convention a' — I. We can now understand why the O 
operators gave correlation functions with sets of poles and were, depending 



on momentum, capable of exciting several special states p5[ . They were, 
in fact, linear combinations of all special operators of a given gravitational 
dimension i.e. spin J. The above construction with the generators (|9.1|) of 



the Woo disentangles the correlation functions. This means that the matrix 



84 



model operators to be identified with the Liouville model special states are 
those defined in ( |9.5|) . 

Let us briefly comment on the meaning of the momentum poles. As 
emphasized in |49[ we should not treat the poles and the | log/i | in the 
correlation functions asymmetrically since the source of the | log/z | is also 
a momentum pole. In fact, all the poles should be thought of as cut off by 
| log/i |. A general N point function (without zero momentum operators) 
would then have | log/i 1^. This proliferation of logarithms was also noted 
in H. 

Let us make some further tests of the correlation functions in ( |9.12| ) and 



(|9.19|) . This will illustrate some very simple consequences of the Woe algebra. 
In fact, the seemingly innocent representation of the puncture and the dilaton 
as /i and t 2 derivatives respectively is a reflection of the W^. Let us give 
a formal argument for this. First the puncture. Write the SU(2) quantum 
numbers of the puncture as J and m, which both will be taken to zero. 
Choose one of the operators in (|9.19|) to be a puncture. We get 



< W X W 2 P >= -(Ji + J 2 



2mi 



2m! — pi ix 



(2/x) 



log/x 



(9.29) 



which, by comparing with the two point function fl9.11p , shows how the 
puncture is represented derivative. The case of the dilaton is equally 

simple. Proceeding as above we find 



< W X W 2 D >-- 



, T T . 2mi 
(J1 + J2)- — + 



4m\ 



Irrtx — pi (2m! — pi) 2 



7T 



;(2/i) 



J1 + J2 



log/i 



(9.30) 



If we introduce explicit t 2 s in the two point function we can write it as 

J1+J2 



< W X W 2 > 



— 2mi 



2/i 



2mi-j?i/(-2t a ) 1 /2 7r (-2t 2 ) 1 /2 \(-2t 2 y/\ 
Taking a (— 2t 2 ) 1 ^ 2 derivative, we arrive at 



log/i I . 
(9.31) 



< W X W 2 D >-- 



( Ji + J 2 + 1) 



2mi 



2mipi 



2mi — p\ (2mi — pi) 2 



7T 



■(2/i) 



J1+J2 



log/i I, 
(9.32) 
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which indeed reproduces ( p.30| ). We must now return to the issue of how 



precisely the dilaton is defined. Recall the original matrix model action 

P J dt[p'\ - ^p 2 - t 2 X 2 } (9.33) 

with (3 dimensionless, £2 having the dimension of energy squared, and p 2 and 
A 2 the dimensions of energy and one over energy respectively. To obtain 
( |9.31|) as a generating functional for dilaton insertions with the above defi- 



P J dt[p\ 



nition of the dilaton we should rescale A and p to make them dimensionless. 
We find 

^H 2 ) 1/2 (P 2 -A 2 )]. (9.34) 

Hence the matrix model dilaton should be represented by (— 2i 2 ) 1 / 2 deriva- 
tives. This is the rescaling eluded to in chapter seven, in the context of the 
Wheeler de Witt equation. 

Let us give some further illustrations in the case of puncture and dilaton 
insertions from the field theoretic point of view. We begin with the puncture. 
Starting with a general correlation function and inserting a puncture does 
not change the Veneziano-like integral which has to be calculated. When we 
insert a puncture, we also must remove one of the screening insertions. The 
only thing which changes is the zero mode part of the calculation. We recall 
the result 

be nuj>-Ae-t _ r(-m)A m . (9.35) 



If we start with r(— m)// m we end up with T(— m + l)^" 1-1 = -mr(-m)/j ra_1 
when we remove a screening insertion. Comparing with fl9.29|) this shows the 
origin of the related factor in the tachyon correlation function. It is a 
consequence of the changing number of puncture screening operators needed. 

For the tachyons the issue of connected or 1PI amplitudes is trivial for our 
case with just one tachyon of differing chirality. There can not be any internal 
punctures just from kinematics. This is no longer the case when we turn to 
the dilaton. The crucial point is that the dilaton can be represented as a t 2 
derivative. Usually, this is precisely equivalent to using the ordinary free field 
contractions, giving Veneziano-like correlation functions. Inserting a dilaton 
in some tachyon correlation function means taking derivatives with respect to 
t 2 (i.e. 1/a'). For dimensional reasons, all tachyon momenta are accompanied 
by a t 2 . Without explicit t 2 s one could write k-j^ for the dilaton. For a dilaton 
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insertion in a nonzero momentum correlation function, the dominating pole 
contribution comes from letting the t 2 derivative act directly on the poles. All 
other terms are clearly less singular. This gives the second term in ( |9.32| ). At 
zero momentum we must also consider the dependence from the t 2 J s which go 
together with the /i's. The latter is a consequence of dealing with 1PI rather 
than connected correlation functions. One could in fact obtain the result 
by considering the explicit combinations of 1PI amplitudes into connected 
ones. To be more precise, if we want to obtain the 1PI amplitude from 
the connected amplitude, we must amputate external puncture legs but also 
subtract of the diagrams with internal puncture propagators. In particular, 
we need to subtract a diagram where a puncture goes off and converts into 
a dilaton. This diagram, therefore, involves a puncture insertion and gives a 
contribution corresponding to a \i derivative. The procedure is illustrated in 
figure 9.3. This is then simulated by an explicit t 2 accompanying the /j's to 
assure the proper subtraction, which corresponds to the first term in ( |9.32| ) 
and only becomes relevant compared to the second term at zero momentum. 
Otherwise we will get one power less of | log/x |'s. We need not restrict 
ourselves to a dilaton among special tachyons. The same reasoning works for 
a dilaton inserted in a general correlation function. 

From these two examples we can conclude that the 1PI nature of the 
correlation functions is very important in the case of zero momentum. 

Finally, it is useful to think of the correlation functions from a group 
theoretic point of view. The symmetry can be used to determine all corre- 
lation functions given the special tachyon correlation functions which may 
be computed by other means. The reason is that all J and m dependence 
of any correlation function is given by some combination of Clebsch-Gordan 
coefficients. For given J's we need the Clebsch-Gordan coefficients of SU(2), 
the 3j symbols, to get the m dependence. To be more explicit let us how- 
ever look at an example, the three point function, to see how the invariance 
properties determine the correlation functions. 

The three point function is obtained by considering coupling (J 1? mi), 
(J 2 , m 2 ) and (J 3 , m 3 ) (with mi + m 2 + m 3 = 0) to (Ji + J 2 + J3 - 2, 0). A 
complication is that there are in general several different channels to sum 
over. This is true already for the three point function. The reason, as we 
have seen, is that we actually should think of the three point function as a 
four point function. The fourth leg carries the excess Liouville momentum, 
i.e. J quantum number, into the vacuum. This is a consequence of the 
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Figure 9.3: Relation between dilaton 1PI and connected correlation func- 
tions. 
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nonconservation of Liouville momentum. Let us use the tachyon three point 
function for normalization. It is given by 



< t iT2 t 3 >= ^ ( ) Jl+j2+Ja -i ! lo L {9 36) 

(2m 2 -p 2 ){2m 3 -p 3 ) tt 

where we have kept the normalization choice of ( |9.21|) . Tachyons 2 and 3 are 
of positive chirality, while tachyon 1 has negative chirality. There are two 
possible channels corresponding to either p 2 = 2m 2 or p 3 = 2m 3 , i.e. 1 and 
2 coming together or 1 and 3 coming together. The group theoretic factor 
in each case is simply proportional to a product of 3j symbols. One for each 
vertex. For the 1-2 channel the factor is 



JM \ ( JJ3J' 
m\m 2 m 3 I \ —m 3 m 3 



(9.37) 



where J = J\ + J 2 — 1 and J' = J\ + J 2 + J 3 — 2. Just retaining the m 
dependence and adding the two channels we find 

(2m 2 - p 2 )(Jim 2 - J 2 mi)m 3 + (2m 3 - p 3 )(Jim 3 - J 3 m 1 )m 2 
= -(2m! - pi)Jim 2 m 3 - (2m 2 - p 2 )J 2 m l m 3 - (2m 3 - p 3 )J 3 m l m 2 (9.38) 



which agrees with (25) after using ( |9.36|) to fix the normalization and J 



dependence. This should come as no surprise since we have just rephrased 
the previous calculation slightly. 

Another convenient way to obtain more general correlation functions is 
through factorization. This is actually already implicit in our previous cal- 
culations. In fact, if we look at ( |9.22| ), we see the complete factorization of 
the tachyon correlation function into a product of three point functions, each 
given by a 3j symbol times a single zero momentum one point function. This 
last piece represents the extra leg in any correlation function which absorbs 
excess Liouville momenta. One may note that these three point functions in 
fact involve states of the wrong dressing. This was also pointed out in [171 



Strictly speaking, the expression in (|9.22|) is just for one channel, the one 



where 1 fuses with 2 then with 3 etc. All channels,however, give identical 
contributions and can not be distinguished. Clearly, the tachyon correlation 
function is consistent with the single Woo factorization result. 
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9.5 Summary 

In this chapter we have exploited the matrix model algebra to derive 
correlation functions. The methods we have used provide a substantial sim- 
plification compared to previous calculations. Also, the correct association 
of special states and elements disentangles the correlation functions in 
a neat way. 

So far we have just looked at the sphere. In the next chapter we will see 
how the picture generalizes to higher genus. 
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Chapter 10 



Higher Genus Correlation 
Functions 



10.1 Introduction 

In the previous chapter we showed how to use the algebra of the c = 1 
matrix model to derive correlation functions of special states. In particular, 
the method is applicable to the special tachyons. In this chapter we will treat 
the higher genus case more explicitly. The motivation is that by deriving the 
results using the W^, rather than just brute force calculations, the crucial 
physics will be more apparent. Hopefully one can in this way obtain useful 
guidelines for constructing a field theory approach to higher genus. In 



similar results as in the previous chapter were obtained independently. The 
authors combine the properties of the system with the powerful path 
integral methods of |58] to calculate correlation functions. This is an useful 



approach for studying the matrix model. We feel however, that the approach 
advocated in previous chapters and used in this chapter, will prove very 
useful for field theory purposes. Furthermore, there are certain subtle issues 
like ordering in the matrix model definition of the special states which can 
not be addressed in the path integral approach. 

An important feature of the matrix model higher genus calculations is a 
deformation of the classical, genus zero algebra with Planck's constant, 
i.e. the genus coupling, as parameter. This is simply a result of replacing 
the Poisson brackets with commutators. 
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As mentioned above, at higher genus we must also face the problem of 
ordering when we define the special state operators. For the special tachyons 
there is no ambiguities, but for the general special state there are. However, 
in the very explicit calculations we will be doing, the question of ordering will 
be important even in the case of tachyons. Not for the final results, but for 
the intermediate steps. We will choose Weyl ordering. Weyl ordering means 
that one takes an average over all possible orderings with equal weights. The 
reasons for this choice are several. It is from many points of view the most 
"natural" ordering prescription. In a path integral approach it would be the 
automatic choice. This is of course connected with another reason, that there 
exist simple formulas for calculations involving Weyl ordering. An important 
question is of course whether this is the ordering prescription which will be 
relevant in the comparison with field theory. It is conceivable that there will 
not exist an unique choice but rather that we have to live with this ambiguity 
which, after all, only amounts to a linear transformation. 

10.2 The Algebra 

We begin our discussion by considering the algebra of the c = 1 matrix 
model. Our conventions throughout this chapter are such that a' — 1. As 
explained in the previous chapter, the theory is organized by a algebra 
generated by the matrix eigenvalue and its conjugate momentum. Classically, 
it is given through the Poisson brackets by 

{W Jl>mi , Wj 2 ,m 2 } = 4i(m 2 Ji - mi J 2 )Wj 1+ j 2 - 1)mi+m2 (10.1) 

where 

W J>m = (A - ip) J+m (\ + ip) J ~ m = 2 J a) J+m a J - m . (10.2) 

For convenience, we have introduced the step up and down operators a) and 
a. As shown in the previous chapter this algebra may be used to organize the 
special state correlation functions in the theory. Parallel work has revealed 
a similar structure in the field theory |5D| , although the precise connection 
between the two approaches is not yet clear. The great advantage of the ma- 
trix model is of course that we may extend the calculation to arbitrary genus. 
The main new ingredient at higher genus is the deformation of the algebra 
mentioned in the introduction. At higher genus, i.e. quantum mechanically, 
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we must also be careful and choose an ordering prescription. As explained 
above we choose Weyl ordering where the new algebra is simply given by the 
Moyal brackets 

1 2 ft d d d d 

[ Wl , W 2 ] = {W U W 2 } M = -sin-(— — - ——)W 1 W 2 (10.3) 
in n 2 op 2 oXi dpi o\ 2 

which allow for simple calculations. The proof is elementary, but we provide 
it for completeness. 

A Weyl ordered operator can be written as 

A(p, q) w = J A{p, q)e- il * Tp+dq) e i{T i> +eil) dpdqdTd6. (10.4) 

If we take the commutator of two such Weyl ordered operators, Aw and Bw, 
we encounter the expression 

^(ti^^ e i{r 2 p+e 2 q)^ = 2 • gin 5( n02 _ T2dl y(Tx+T2)p+i(6i+e 2 )q) _ ( 1Q 5 ) 

The sin is then rewritten as 2zsin|(^-^- — ^gfj) by partial integration. 
The T\ — t 2 and Q\ — 9 2 integrations give 5 functions which fixes pi = p 2 and 
Qi — <l2 respectively. From this the result follows. 

An equivalent way of exhibiting the quantum deformation is to focus on 
( |10.5| ), which in fact is a very old result, p2| , and inspired by this introduce 



a generalized loop operator [jrj| . Instead of (|8.6| ) we define 

w(k, I) = e kp+lx , (10.6) 
where p and A are the conjugate variables. We then have 

[w(ki, h),w(k 2 , l 2 )\ = 2i sin ^(hl 2 - l 1 k 1 )w(ki + k 2 , h + l 2 ) (10.7) 

with h — > giving back a W^. Interestingly, it can be shown [|j] that ( |10.7| ) is 
a representation of SU(N) with h = 1/N. This is reminiscent of the original 
unitary symmetry of the matrix model. 

Let us give a simple example of a higher genus commutator. Using the 
Moyal bracket it immediately follows that, for h — 1, 

oo i 2g 
m _fni _ \ ^ i TT/Vn _ — h\ ( n m-2g-l Jn-2g-l 



\a , a 



£ 9W9 IK™ " *)(» ~ k)(a m - 2 ^a^ 2 ^) w . (10.8) 
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The continuation to the upside down harmonic oscillator gives an extra factor 
of i and an extra (— l) 9 . This result will be extensively used in a later section 
where we calculate tachyon N-point functions. 

For clarity let us calculate explicitly a two point function to all genus. 
We choose the correlator between spin J = 3/2, m = 1/2 and J = 3/2, 
m — —1/2. To do that we need to calculate < W 2 $P >■ This is easy. We 
have 

< wvp >= --imf; ^ 2 ; 1 . (io.9) 

K n=0 h n - t 

The extra term —1 comes from Weyl ordering and will be thoroughly dis- 
cussed in the next section. If we keep only terms nonanalytic in to this 
reduces to 

< W 2fi P >= (4*g - 1) < PP > . (10.10) 

To evaluate our two point we use the Moyal bracket continued to the upside 
down harmonic oscillator to calculate 

[^3/2,1/2, W 3 /2,-l/2] = 2(6^2,0 + 4ty ,0) (10.11) 

( plOD and (|TaT0| ) then finally give 



< ^3/2,1/2^3/2,-1/2^ >= (24tjj — 2) < PP > j—. (10.12) 

The same procedure may be used to calculate arbitrary correlation functions. 

In the previous chapter, it was shown how to calculate special state cor- 
relation functions using this algebra. The number of operators in the cor- 
relation functions are reduced one by one using the algebra. The last step, 
however, involves the calculation of a remaining zero momentum one point 
function. This is where the cosmological screening charges sit. Let us there- 
fore turn to the issue of calculating such objects. 



10.3 The One Point Functions 

To obtain the zero momentum special state one point functions at higher 
genus, we need to be careful. From the algebra we understand that they are 
given in terms of (H n )\y- These are just powers of the Hamiltonian. But it 
is important to remember the Weyl ordering. This will make the one point 
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functions quite complicated. We need to rewrite the Weyl ordered powers in 
terms of something simpler. The objects with simple correlation functions 
are clearly the H n , which give 

t%<PP> (10.13) 

since the resulting in the trace can be replaced by up to terms analytic 
in to. To obtain the {H n )w one point functions we therefore need to perform 
the Weyl ordering and express the result as a polynomial in H . One way to 
do this is to use the standard formula 



m 



(a m a^ n ) w = 2~ m J2 ( ? a? m - l a n (J l . (10.14) 



1=0 



See e.g. [Q. If we put m = n and take the expectation value in the eigenstate 
r we get 

n / \ n— 1 

2 ~ n E / )U(r + l-k), (10.15) 

1=0 \ / k=0 

and finally 

n / \ re— 1 i 

(H n ) w = 2""E " H(H-- + l-k). (10.16) 
i=o \ 1 J k=0 / 

This is unfortunately a very complicated polynomial. In the next section we 
will explicitly calculate some tachyon correlation functions. To do so we need 
explicit expression for the coefficients in (|10.16|) . By some tedious algebra, 
which is greatly simplified if Mathematica is used, the following expression 
can be obtained 

(H n ) w = H n + ±(2n - l)n{n - l)H n ~ 2 

+^(20n 2 - A8n + 7) nLo(" ~ (10.17) 

+ 29oko( 280n3 - 1596n2 + 187An - 93 ) n!= (« - k)H n ~ 6 + ... 

The continued case is obtained by sending H n ~ 2k — > (—l) h H n ~ 2k . This 
is enough up to genus 3 for the general correlation function. In appendix 
10A, the expansion is shown up to genus 7. These expressions are really 
remarkably complicated, and responsible for the complicated structure of the 
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higher genus correlation functions. The piece coming from the algebra 
is much simpler. 

For clarity, let us give an example. We choose J = 4. We find 



where we note resonance contributions up to genus 2, but no higher. The 
same general structure is clearly found for the tachyon correlation functions. 
Poles are present only up to some genus. The precise value depends on the 
scaling (i.e. momentum). 



We now need to combine these results with those given by the deformed 
algebra to get the tachyon correlation functions. This is, in principle, rather 
simple. In fact, one can write down a set of diagrammatic rules which give 
the correlation functions. 

By the variational procedure introduced in the last chapter one can show 
that in the case of (N — 1,1) amplitudes it is possible to arrange things in 
such a way that only one basic graph is needed. There exist other choices 
too, involving sums over several graphs, but this one is clearly the simplest. 
The graph, shown in figure 9.2, has a T~ at one end and the necessary one 
point function at the other. In between are the T + distributed in some order. 
The order does not matter, but we should pick just one ordering, no sum. 
To see this, one can start with the two point function of a T + and a T~. 
All subsequent T + variations will only give a nonzero value when they hit 
the — leg. If we are considering a genus g correlation function, we must 
partition the handles in all possible ways among the vertices and the one 
point function. Each such diagram is then given a certain factor according 
to the following rules. 

For each vertex write down a factor f 9k {pk) given by 





(10.18) 



10.4 Tachyon Correlation Functions 



fgM 




(Pk ~ l)(Pk - 2)...(p fc - 2g k ), 



(10.19) 



2^(2g k + l)\ 
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where pk is the momentum of the fc'th (positive) tachyon and gu the genus 
of the corresponding vertex. This is obtained from ( 10.8Q . Strictly speaking 
the formulae are only valid when p = 2 J, i.e. for special tachyons. The one 
point function is given a factor h g according to 



ho(p) = 1 

h 1 (p) = ^(2p-l) 

h 2 (p) = ^(20p 2 - 48p + 7) 

= 2903k ( 28 °P 3 " 1596 P 2 + 1874 P - 93), 



;io.2o) 



where the argument p is 



N 



\Pk I -N + 2(~g-g) + l=\ Pl \ -N + 2(g-g) + l. (10.21) 



k=l 



g is the genus associated with the one point function. Clearly, J2k=2 9k+9 = 9- 
Some additional h g polynomials may be read off from appendix 10A while 
the general case is generated by (|10.16|) . In both (|10.19|) and ( |10.20| ) we have 
extracted a common factor depending only on the momentum of the single 
negative chirality tachyon. 

Using this technique, the two point function would be calculated as 

p 2 fo(p)h (p - l)- 2 — + p 2 (p - l)(p - 2) 



p 



t: 



p-3 



x [fi(p)h Q (p - 3) + / (p)/ii(p - 1)] -5— + p\p - l)(p - 2){p - 3)(p - 4) 

P — 6 



x [f2(p)h (p - 5) + h(p)h x (p - 3) + / (p)/i2(p - 1)] 



p— 5 



+ ... (10.22) 



p — 5 

up to genus 2. In particular, if we concentrate on genus 2 the relevant 
calculation involves 



p 2 {p- l)...(p-4) 



1920 

^-l)(p-2)(2,-7) + ^ 



(p-l)...(p-4)+ 
(20 (p- l) 2 -48(p 



7) 



p 2 (p - l)...(p - 4)— -(3p 4 - 10p 3 - 5p 2 + 12p + 7). 
5760 



10.23) 
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Although the algebra is somewhat tedious, the generalization to higher 
genus and N-point functions is straightforward. As an example the three 
point function at genus 3 is computed in appendix 10B. Here, we will simply 
state the results for N-point functions up to genus 3. 

P)N-3 i N f)N-l 

+P1- 1 __fr„2_ ))i _ i \fi ^1-1 



dt N-3"o 24 V ^f* yi J dt%- 10 

N N N N N 



+^(3E^+io E p^ 2 -iopiE^ 2 -5( E (ft-p,-) 2 -(Jv-3)Erf) 

i=2 i,j=2,i>j i=2 i,j=2,i>j i=2 



QN+1 

-i N N N 

' (9E^ + 63 E P?P} + 210 E P?^Pi 



2903040 i=2 ij=2,i>j i,j,i=2,i>j>i 

N N N N 

-Pi(63Ep' + 210 E P?P?)+42EPi+210 E 

i=2 i,j=2,i>j i=2 i,j=2,i>j 

N N N N 

+210 ]T P^p 2 +Pi(217^ Pl 2 -70 £ ftp,)- 294 £ 

i,j,l,i>j i=2 i,j=2,i>j i,j=2,i>j 

- 205p 1 -93)^ rF 3^- 1 . (10.24) 

This is to be multiplied with the appropriate factorized pole factor. The 
result may be compared to other calculations of the correlation functions. 
Clearly this is a very complicated expression which practical application is 
probably very limited. The main point is, however, not this expression in 
itself, but rather the way in which we have derived it. 



10.5 Summary 

We have seen how to calculate special tachyon correlation functions in the 
matrix model using the deformed algebra. The algebra nicely organizes 
the correlation functions even at higher genus. 

A natural challenge to the field theory approach is to reproduce this de- 
formation of the algebra, which we have described and found to be important 
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for the structure of the correlation functions. One would expect that the de- 
formation is related to contour integrals of the currents which when pulled 
back over the surface, away from any insertions, fail to yield zero since they 
get caught up around handles and cannot be contracted to points. 
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Appendix 10A 

Weyl ordering of powers of the Hamiltonian for computations up to genus 7: 
(H n ) w = H n + ±(2n - l)n(n - l)H n - 2 



n-4 



+ ^(20n 2 - A8n + 7) Ul =0 (n - k)H 

+ 2903040 ( 280 ™ 3 - 1596n2 + 18Un - 93 ) nLo(« - k)H n ~ 6 

^2Qo(2800n 4 - 29120n 3 + 85544n 2 - 68368n + 1143) I\l=o(n - k)H n ~ 



1393459200 



+ 367873228800 ( 12320n5 ~ 203280n 4 + 1100528n 3 - 2255352n 2 
+ 1422434n - 7665) IlLo( n - k)H n - w 

+ 24io30539 1 509760oo (H211200n 6 - 269068800n 5 + 2312950640n 4 - 8743967232n 3 
+ 14194825268n 2 - 7658714592n - 12730293) HHoi™ ~ k)H n - 12 
+ 578473294823424000 ( 32 3 2 0n 7 - 105385280n 6 + 1315394080n 5 - 7873930064n 4 
+23465896792n 3 - 32333229596n 2 + 15636417798n - 7740495) Uli (n - k)H r - 



100 



Appendix 10B 

In this appendix we calculate the tachyon three point function at genus 3. 
As explained in the main text we should write down all possible partitions 
of the three handles among the vertices and the one point. This gives 

[/s(Ps)/o(P2) + /o(P3)/s(P2) + /2(ps)/l(P2) + fl (P3 ) h (P2 ) ] &o(pi ~ 8) 
+ [fMhip-l) + /o(P3)/ 2 (P2) + /l(P3)/l(p 2 )] ^l(Pl - 6) 

+ [/i(p 3 )/o(P2) + /o(p 3 )/i(p 2 )] h 2 (pi - 4) + /o(p 3 )/o(P2)/i 3 (Pi - 2). (10.25) 
Using the expressions (|10.19|) and (|10.20| ) we find 

1 



2 6 7! 



[(p 3 _l)...(p 3 _6) + (p 2 -l)...(p 2 -6)] 



+^^y [(ps - l)-(ps - 4)(pa - l)(p 2 - 2) + (p 2 - l)...(p 2 - 4)(p 3 - l)(p 3 - 2)] 

+ ( 2^5! [(P3 " 1) - (P3 ~ 4) + (P2 " 1} - (P2 " 4)] 
l(Ps - l)(ps - 2)(p 2 - 1)(P2 - 2)] - 13) 

+ 2^! [(P3 " 1)(P3 ~ 2) + (P2 " 1)(P2 " 2)] 57^0 [ 2 ° (Pl " 4)2 " 48(Pl ~ 4) + ? 

+ ^5^( 28 °(^ " 2 ) 3 " 1596 (^ - 2 ) 2 + 1874 ^ " 2 ) " 93 )> ( 10 - 26 ) 
which is nothing else than, if we use momentum conservation, 

1 



:(9(pE + Pi) + 63(p 2 2 p| + p^) - Pl (63(p* + p|) + 210p*p|) 



2903040 

+42(p* + p£) + 210(p 2 p^ + pin) + pi(217(p^ + pi) - 70p 2 p 3 ) 

- 294p 2 p 3 - 205pi -93), (10.27) 



which indeed is a special case of (|10.24j) 
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Chapter 11 

Comparison with Liouville 
Theory 



11.1 Introduction 

In this chapter we will make some simple comparisons between the matrix 
model and field theory calculations. The first and most obvious observation 
is of course the presence of the special states in both theories. This is in 
fact one of the main results of the previous chapters. Furthermore, there is 
a Woo-algebra organizing the correlation functions in both cases. All of this 
lend support to the accepted view that the theories are, in fact, identical. A 
very nice and explicit connection at the level of the algebra was given in [ITS 



Clearly, one would also like to have a precise comparison of the correlation 
functions. As far as tachyons go, we have seen the agreement to be perfect. 
At least on the sphere, the only case for which field theory calculations in 
general exist. When comparing results for the special states, additional care 
must be taken since we are automatically sitting on the discrete momentum 
poles. This can be especially tricky in the field theory approach. 

An important issue to decide is whether the correlation functions one is 
calculating are correlation functions of states 

W(z)W(z), (11.1) 

or currents 

W(z)0(z) + O(z)W(z). (11.2) 
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The O's are ghost number zero spin zero fields needed for well defined closed 
string currents. They build up the "ground ring". For a discussion see 
[ |7§| , |S0j. As for (A ± p) 2J there cannot be any disagreement. Clearly, they 



give correlation functions of tachyon states. After all, as explained in [|78| 



there are no currents corresponding to the special tachyons! The needed 
elements O of the ground ring do not exist. With this identification for the 
tachyons it would be very unnatural to associate anything else than special 
states to the other correlation functions of chapter 8. 

Let us try to be more explicit. We will first look at a very simple example 
involving the dilaton, and then proceed to a more general case. 



11.2 A Dilaton Comparison 

There are some very simple examples of correlation functions easily com- 
putable just using Liouville notions and no matrix model techniques. These 
are correlation functions involving the dilaton. We will in fact be able to ob- 
tain some results to all genus simply from dimensional arguments. Consider 
the Liouville partition function (or space time free energy) 

E(A) = lim - [ PIP^-s/(- 2WOX9 " x+ i ao **^W A ^, (11 .3) 

R-*oo R J 

where t 2 = — ^7, Q = 2V2, a = — Vz and R is the radius of the target 
space for the matter field X. A is the world sheet cosmological constant, 
dimensionless from the point of view of space time. The only dimensionful 
quantities are R and a'. In the noncompact case we have in fact only a' at 
our disposal. From dimensional grounds and KPZ scaling we must have 

E(A) g ~ (-t 2 ) 1/2 A 2{1 - 9) (11.4) 

at genus g. E(A) is the generator of connected amplitudes (in space time). 
Let us perform a Legendre transform, as we did in chapter seven, to obtain 
a generating functional for 1PI amplitudes with respect to the puncture, i.e. 
the zero momentum tachyon. This means taking away any pinches. We have 

£7(A) = A//-rM, (11.5) 

with 

/^=|f = H 2 ) V2 (A + ...)• (H.6) 
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\i has the dimensions of energy. Hence 

r(//) ~ (-* 2 )*< 2 »-V~ 2fl (n.7) 

In the 1PI generator T we should of course regard // as independent of t 2 . 
Since t 2 derivatives should generate dilaton insertions we find the following 
1PI amplitude relation 

1 n 

< 2 ...0 2 > 9 = —— U(2g + 1 - 2p) <>„ (11.8) 

l^2j p=1 

which is identical to what was obtained in chapter seven, using the matrix 
model recursion relations (generalizations of the zero momentum Wheeler de 
Witt equation). As noted there the dilaton one point function involves a 
factor 2g — l rather than the expected 2g — 2. From the above it is clear that 
this discrepancy is simply due to including the overall (— £ 2 ) 1,/2 i n ( 11.4Q . 



In general we cannot get away with a simple trick like this, more work is 
needed! 



11.3 More General Correlation Functions 

To obtain the general special operator correlation function in the Liouville 
theory one would like to use the group theoretic information provided by the 
Woo or, for given spin J, the SU(2) symmetry. The states in the Liouville 
theory are given by combinations W(z)W(z) of the Liouville theory version 
of the special states W(z). Given this it is tempting to believe that we have 
a representation of a W^ x W^ symmetry (left times right). However, this 
is in general not correct. In the uncompactified case the left and the right 
moving states must be the same. The symmetry group is broken down to just 
the diagonal subgroup. This is achieved in two steps. First the gravitational 
dressing must be the same for left and right, otherwise we would be unable to 
screen using the cosmological constant which treats left and right in the same 
way. This means that we always must have the same spin J for left and right. 
We get a reduction to the diagonal of the piece transverse to SU{2) x SU (2). 
This is true even for the compactified case. If we in addition are considering 
the uncompactified case, the left and right moving momenta must be the 
same and therefore also the m quantum numbers. Consequently, we just 
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have a representation of the diagonal Woo- This is in precise agreement with 
the matrix model, where we indeed only see one Woo- 
There is however an apparent paradox here. If we would use the free field 
contractions in computing the correlation functions the results would seem 
to disagree since from this point of view left and right are still independent. 
From such a calculation you would expect to get a different result, all Clebsch- 
Gordan coefficients squared, one from the left and one from the right. We 
will however show that the results in the end turn out to be consistent. 
Let us begin by considering the tachyon correlation function as computed 
2M. As we have seen the result is in complete agreement with the matrix 



m 



model results. On the other hand we have seen how the matrix model orga- 
nizes its correlation functions using a single Woo- Let us consider the Liouville 
calculation more carefully. The result of p4j is obtained through arguments 



of analyticity and symmetry. In particular the by now well known factor- 
ized product of gamma functions is found PjL |36|j with a certain unknown 
coefficient independent of the particular momenta. This coefficient is then 
determined by sending all the momenta, except three, to zero. This reduces 
the expression to a three point function with N — 3 extra punctures. Since 
the three point function is possible to evaluate directly, the general result 
follows. The extra N — 3 punctures is simply represented as n derivatives. 
This is the Liouville derivation of the expression ( |9.21| ). The important point 
is that the use of a \i derivative for inserting a puncture is a consequence of 
having just one Woo-, see ( |9.29|) ! This means that the calculation in [f24 



automatically incorporates this feature. 

For the more general case with nontachyonic special states, we return to 
the recursion relation (|9.20| ). Let us redefine the fields according to 



AT 

Wj, m = Wj, m . (11.9) 

2m — p 

The recursion relation then takes the form 

<f W TTT >- (Ji-mi)(J+Ji-l) c & £ f > 

< l.J,JVVJ 1 ,m 1 11 i- Ji,Ji >— j < VV j +Jl - ly j +mi 11 lj uJi > 

(11.10) 



In |47] these very same recursion relations were obtained in the case J = m = 
1/2 using Liouville methods. The coefficient in front of the right hand side 
was shown to be of the form (2 J\ — 1)C 2 , where C stands for the appropriate 
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Clebsch-Gordan coefficient. The first factor comes from comparing with the 
purely tachyonic case where the answer is obtained from a simple Veneziano 
like integral. To see the agreement one uses the Clebsch-Gordan coefficients 



of the special operator algebra as obtained in [48 



yj 3 ,m 3 _ A(J 3 , m 3 ) 

A(J 1 ,mi)A(J 2 ,m 2 ) 



(±:::u,,,,, = — T „ - - j^o, (n.n) 



where 

A(J,m) = -- [(2J)!(J + m)!(J-m)!] 1/2 . (11.12) 



At J = m = 1/2 one finds C 2 = &~p> which then leads to (pTT0|) . This 



is an important check on the equivalence between the Liouville and matrix 
model approaches. An everywhere present difficulty in these comparisons is, 
however, the fact that we are really sitting right on the momentum poles. 
Clearly one needs to carefully regularize all expressions. 



11.4 Summary 

We have investigated the structure of special operator correlation functions 
in c = 1 quantum gravity. Due to the presence of a symmetry the 
calculations become very simple. We have also investigated the connection 
between the Liouville and the matrix model, indicating the agreement for 
the correlation functions. 

An important point is the existence in the uncompactified Liouville as 
well as matrix model formulation of c=l of just one Woo. From the Liouville 
point of view this is somewhat obscure since the operator product expansion 
and its Clebsch-Gordan coefficients seem to give a structure corresponding 
to two Wqo's. Fortunately, the final outcome of the explicit calculations are 
identical. Further work is however needed to establish the full equivalence. 

A part of the problem is that many of the calculations are so ill defined. 
The reason is that we are sitting right on the discrete momentum poles. 
Especially in the Liouville theory this is a big technical problem. Often we 
must rely on guesswork concerning ill defined analytical continuations. It 
is very doubtful if many of the results would have been obtained correctly 
without knowing the answers in advance, given by the much more powerful 
matrix model. 
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Chapter 12 
Conclusions 



In this thesis we have studied the two dimensional noncritical string. It is a 
remarkable example of a string theory where we are capable of calculating 
everything. This is, of course, thanks to the matrix models. 

The study of the two dimensional string theory has so far been of a quite 
technical nature. The matrix model has been used at length to extract all 
possible correlation functions and to study the spectrum. One of the main 
results of this thesis is in fact a clear verification of the presence of the special 
states in the matrix model. 

Clearly, this should not be a goal in itself. After all, the reason why we 
are studying this model in the first place, is that we want to learn something 
about string theory and quantum gravity. There are several different ways 
to do that. So far, the most important contribution of the matrix model 
has been rather indirect. It has been to act as a source of inspiration for 



new developments in field theory. As is clear from many recent works, [46 
|5T| [771 , Hl> H> there has been considerable progress in the understanding 
of the field theory techniques. Perhaps, eventually, this will help us in the 
study of higher dimensional theories. A somewhat rebellious mind would 
perhaps argue that it is instead time to abandon the old fashioned customs 
and stick with the matrix model. Unfortunately, this is clearly not possible 
since the matrix models seem to have a rather limited range of applicability. 
In particular there is the c = 1 limit. It is however an important lesson to be 
learned, that there, in some cases, can exist totally novel techniques which 
might lead to exciting breakthroughs. 

On the other hand, in the study of the specific two dimensional string 
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theory, remarkably little has been done with respect to the understanding of 
the target space picture. Not only is the relation between the Liouville mode 
and the matrix model eigenvalue rather obscure (recall the relation through 
the loop length in the Wheeler de Witt equation) we have in addition very 
little control over the target space background. In the field theory approach 
we can choose to put in a black hole etc. We would like to be able to exercise 
a similar control in the matrix model. Only after we have learnt how to 
do this, we can truly use the matrix model to learn about stringy quantum 
gravity. 

There are also several issues regarding the back reaction of the background 
geometry against the loop corrections to the vacuum energy. The naive loop 
corrected f3 functions clearly indicate that there should be a highly nontrivial 
influence. From this point of view it is not clear how a flat space time 
can remain a solution at loop level. However, arguments in J73J seem to 
indicate that the two dimensional background geometry is quite rigid and 
not influenced in this way, but rather frozen. If this is true, there is a risk 
that the things we might learn which are relevant for higher dimensions, 
where this certainly is not the case, are rather limited. 



Other work point in a perhaps more promising direction, [H, [TO, [TJ, [70 



76fl . There, a one loop corrected string inspired field theory is studied for 



black hole solutions with subsequent Hawking evaporation. An interesting 
elaboration would clearly be to assume an underlying string theory and to 
study its consequences. A better understanding of the matrix model could 
be of invaluable help in this respect. 

Clearly, there are still a lot of things which remain to be done in the matrix 
model approach to the two dimensional string. We need to take advantage of 
its simplicity and learn everything we can. But it is also important to look 
forward, to more physical models. Although the matrix model techniques as 
they are now seem to be insufficient for such a task, there could clearly exist 
generalizations which might help us in our quest. 

The most important lesson we might learn from the matrix models is the 
importance of keeping an open mind. The field theory formulation might not 
be the only way of dealing with strings or whatever it is that we feel deserves 
our attention. Finally, however, it is important to keep the remarks made in 
the introduction in mind. We should not become too carried away with the 
progress made. Even if our mathematical capacity is continuously increasing, 
this does not mean that our understanding of physics is any deeper. 
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String theory is an exciting possibility, but whether it is the right answer 
only nature itself can tell. 
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